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ABSTRACT 


Title of Dissertation: Rapid X-Ray Variability of Active Galaxies 

Allyn F. Tennant, Doctor of Philosophy, 1983 

Dissertation directed by: Ellhu A. Boldt 

Adjunct Professor 
Department of Physics 


Active galactic nuclei are luminous sources of X-rays. It has 
generally been assumed that the X-rays are generated within 10 
gravitational radii from the central object. Research in this thesis tests 
this idea by making a very sensitive search for rapid (< 1 day) X-ray 
variability from active galaxies. 

To perform this search one has to develop statistical techniques to 
separate true source variability from noise. Methods, tested with Monte 
Carlo simulations of the data, are presented which allow this separation 
when the data bins have unequal length and which determine the time scale 
of source variability given an observed source variance. 

The in orbit performance of the detector is quite good. A small 
signal caused by particle contamination is greatly reduced if one only uses 
data taken at low values of Mcllwain L. A study of about 50 observations 
of ’blank sky' shows another signal at a level of 1% of the total 
background. This signal is due to different parts of the X-ray background 


being sampled due to spacecraft jitter and so measures the spatial variance 
of the sky. A 'serendipitous* burst, discovered during a blank sky 
observation, appears to be a flare from a nearby galactic source. 

Observations of 38 different active galaxies show no evidence for 
rapid variability for the vast majority of the objects. Three objects 
which do vary show a time scale consistent with one day. Only the 
observation of NGC 6814 shows a time scale shorter than one day, and this 
object shows factor-of-two changes for periods as short as 2 minutes. 

The non-vari ability of most objects Indicates that the X-ray 
producing plasma Is either stable or large. A nonthermal relativistic 
electron population would explain the X-ray spectra. Since these electrons 
cannot be gravitationally bound, the X-ray plasma will fill a large volume 
relative to the size of central object. Occasional flares are produced 
wher, a new burst of particles is Injected Into the surrounding cloud. In 
this picture NGC 6814 would be dominated by the variable emission from the 
Injection mechanism instead of the more constant emission from the 
surrounding cloud. 
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1. INTRODUCTION 

1.1 Chapter Overview 

In this chapter I provide a general Introduction to the properties 
of active galactic nuclei. I begin with a general description of the 
Important observed properties for the non-expert. Optical observations 
Indicate a strong source of Ionizing photons (hv > 13.6 eV). It js quite 
likely that the observed X-rays are related to this Ionizing source. The 
X-rays are generally assumed to be produced very near to the "central black 
hole". This Idea Is based more on belief than on fact and so needs to be 
tested. The research presented In this thesis makes such a test by 
performing a very sensitive search for rapid {< 1 day) variability. 
Observations of rapid variability could Indicate a compact source of X-rays 
whereas the lack of variability could Indicate extended emission. The 
chapter ends with a discussion of other observations of rapid variability 
both In the X-ray band and at other wavelengths. 

1.2 Properties of Active Galaxies 

Seyfert (1943) pointed out that although many galactic nuclei 
contain emission lines, there is a small class of objects which show strong 
Balmer lines In addition to lines from highly excited states. Over the 
years galaxies showing such lines have come to be called Seyfert 
galaxies. The strongest optical lines are from the Balmer series, twice 
Ionized oxygen [OIII] at 5007 A and once ionized nitrogen [Nil] at 
6583 A. A wide range of ionization Is observed ranging from [Old, [NI] and 
ESII] to [Ne V], [Fe VII] and [Fa X] (Osterbrock 1981). These forbidden 
lines typically have widths (FWHM) ranging from 500 km/sec to 1000 km/sec 
(Khachlklan and Weedman 1974). These lines are called the narrow lines 


since they are often quite narrow when compared to the hydrogen 
recombination lines. 
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The Balmer lines are very strong and can have widths ranging up to 
5000 km/sec. In 1971 Khachlklan and Weedrmn proposed dividing Seyfert 
galaxies Into two classes based on the widths of the permitted lines. In 
Seyfert 2 galaxies the permitted lines have roughly the same width as the 
forbidden lines. In Seyfert l's the permitted lines have widths ranging 
from 1000 <m/sec to 5000 km/sec. 

Seyfert galaxies also have very bright nuclei. As Weedman pointed 
out In his 1977 review of Seyfert galaxies, galaxies with strong broad 
emission lines always have bright nuclei. Some galaxies with bright nuclei 
do not show Seyfert characteristics. However, Weedman points out that If 
the nucleus Is bright relative to the disk, then It Is very likely to be a 
Seyfert. There is a wide range of nuclear luminosities. Yee (1980) 
published the non-thermal luminosities of a large number of active 
galaxies. Although a large part of his non-thermal emission could be due 
to lines, it does represent a good estimate for the luminosity of the 
nucleus. These luminosities ranged from 3 x lO^ 1 ergs/sec up to 3 x 10^ 3 
ergs/sec for the Seyfert 2's. For Seyfert 1 galaxies the range was 3 x 
10 41 to 3 x 10^ erg/sec. The sun emits 2 x 10 33 ergs/sec (Integrated) and 
the luminosity of a bright galaxy equals 10 11 suns. Therefore, a large 
fraction of the total light emitted by an active galaxy can come from the 
nucleus. 

Radio astronomers discovered that some galaxies can be strong 
sources of radio emission. Radio emission from <_ ' xles generally comes 
from three components, with two components forming ■ d e lobes on opposite 
sides of the galaxy. These lobes can be quite large and are bright at low 
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frequencies (MHz). For Cen A, a bright nearby radio galaxy, the radio 
lobes are about 1 Mpc In length and the radio spectrum Is a power law with 
an Index that ranges from 0.9 to 0.6 over a frequency range of 19.3 Mz to 
500 Mz (Cooper, Price and Cole, 1965). The nucleus of a radio galaxy Is 
also observed to emit radio waves at high frequencies. For this work we 
will only be concerned with the nucleus. However, observations made at the 
Very Large Array (VLA) have shown that In a radio galaxy the nucleus Is 
often connected with the external lobes via jet like structures. This 
provides clear evidence that the nucleus also powers the looes. 

Optical spectra of the nuclei of radio galaxies show that they come 
in two types: those with broad lines and those with narrow lines. Thus the 
spectra of radio galaxies greatly resemble the spectra of Seyfert galaxies, 
although In a statistical sense It Is possible to distinguish between the 
two classes (Osterbrock, 1978). The properties of the associated galaxies 
are quite different. Radio galaxies almost always occur in elliptical 
galaxies whereas Seyfert galaxies tend to be spirals. Wilson et al . (1980) 
and Wilson and Willis (1980) observed that some Seyfert galaxies do contain 
very small radio lobes. Wilson and Willis argued that the galaxies 
containing lots of gas and dust (l.e. spirals) do not produce giant radio 
lobes since the gas and dust blocks the energetic particles from reaching 
Intergalactic space. 

By the early 1960's radio astronomers had discovered a large number 
of radio galaxies. Some radio sources did not appear to be associated with 
galaxies but rather with star-like objects having rather unusual optical 
spectra. In 1963 Schmidt Identified the lines in the spectra of 3C 273 as 
being typical lines from an active nucleus but at a redshlft z 
° AX A * .16. Shortly thereafter Greenstein and Matthews (1963) identified 
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the redshlft of 3C 48 to be at a z = .37. In an expanding universe 
redshlfts correspond to distance and so quasars are generally assumed to be 
very distant active galaxies. Over the years many quasars have been 
discovered and a couple have been found with z > 3.5. In the 1960's It was 
discovered that quasars can be either radio loud or radio quiet. The 
recombination lines observed In quasar spectra can be either broad or 
narrow In the same way that Seyfert spectra can have either broad or narrow 
lines. 

1.3 The Standard Model 

It Is useful to have In mind the general picture or model of an 
active nucleus, as this allows greater understanding as to where various 
observations fit Into the total picture. It Is typical to divide the 
nuclear region Into three parts: 

1. There Is a "point source" at the center which Is the power 
source for the nucleus. This source, sometimes called the central engine, 
radiates large quantities of ionizing radiation. It Is generally assumed 
that the ultimate source of power Is gravitational energy release. 

2. Surrounding the central engine and extending out to several 
light days for the low luminosity objects and several light years for the 
high luminosity objects, Is the broad line region. This region consists of 
dense (10 9 - 10 10 cm" 3 ) fast moving clouds at a temperature of ~ 10 4 K. 
Photoionization of these clouds produces the broad lines and colllslonal 
deexcitation suppresses emission from the forbidden lines. 

3. External to the broad line region Is the narrow line region 
which can extend to a radius of several thousand light years. The clouds 
In this region are less dense ( 10 2 - 10 6 cm" 3 ) and so can reradiate 
forbidden lines when photolonlzed. 
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Thus In the standard picture all active nuclei are very similar. 
Narrow line objects are merely active galaxies without broad line clouds. 
Quasars are distant and hence luminous forms of active galaxies. Radio 
galaxies are active galaxies that allow the relativistic particles to 
escape. And finally, BL Lacs may be radio galaxies In which the escaping 
particles are coming roughly directly towards the Earth. 

1.4 Timing 

In the standard model the lines are produced as a secondary source 
of radiation. Although emission lines provide a very good diagnostic of 
the environment In which they are produced, this environment turns out to 
be the broad and narrow line regions. Hence, the lines tell us relatively 
little about the central engine. In order to stucty the source of the power 
one needs to consider the continuum. 

In the standard picture it Is popular to assume that X-rays are 
produced deep in the potential well near a supermasslve object (see Rees, 
Begelman and Blandford 19ET for a recent discussion). If the X-ray flux 
from this dense, gravitationally confined plasma Is observed to vary, then 
the shortest time scale will be on the order of the light travel time 
across the innermost stable orbit in a Schwarzschlld geometry. This is 
given by 


Ax ~ R/c ~ 6 GM/c^ ~ 50 M g sec 

where M g is the mass of the central object in terms of 10 6 solar masses (M 0 
= 1 solar mass = 2 x 10 33 gm). Since these objects are not expected to 
exceed the Eddington limit, we expect Mg to He In the range of .1 to 1000 
for luminosities ranging from 10 43 to 10 47 erg/sec. Thus the relevant time 
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scales range from ~ 10 sec to ~ 1/2 day. Lack of variability on these time 
scales could be an indication that the X-ray emitting plasma is not 
gravitationally contained. 

Fabian and Rees (1979) (see also Cavallo and Rees 1978) have 
proposed a general constraint on the shortest time scale for a given change 
in luminosity. They assumed that a spherically symmetric, homogeneous 
cloud of stationary matter is converted Into radiation with an efficiency 
of n» The shortest observed time for the change to take place Is given 
when the optical depth of the cloud is unity. This gives the minimum 
time ~ 5 L^ 3 /q sec where L43 Is the luminosity In units of 10^ 
ergs/sec, 

1.5 Other Observations of Rapid X-Ray Variability 

An early report or variability on a short time scale was made by 
Winkler and White (1975). They reported that the X-ray flux of Cen A 
Increased by a factor 1.6 over a 6 day observation, which is considered to 
be a mild change by "modern" standards. 

Delvallle, Epstein and Schnopper (1978) found evidence for a 25% 
step Increase In flux from Cen A. They reported a 2 sigma uncertainty In 
the length of the step of 0 to 5 hours. Their observed count rate was 
quite low and so the error for a single time bin was about 25%. This 
result could have been affected by systematics such as nearby sources (see 
Marshall and Clark 1981). However, Lawrence, Pye and Elvis (1977) have 
shown that Cen A Is continuously variable over extended periods of time 
with a time scale of about 1 day. 

The second brightest active galaxy, NGC 4151, also has a long 
history of X-ray variability. Elvis (1976) reported that NGC 4151 flared 
by a factor of 1.7 in less than 3 days. Tananbaum et al. (1978) reported a 
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factor of 10 change on a time scale as short as 730 sec. I discuss thin 
observation In Appendix F and point out that the 700 sec result Is probably 
In error. Mushotzky, Holt and Serlemltsos (1978) reported a flare with a 
factor of two change observed over 1.5 days. Lawrence (1980) finds such 
flare-like events are rather common In Ariel 5 data extending over several 
years. Although some of the observed flux variations on longer time scales 
are caused by variations in the absorption column (Barr et al. 19775, there 
are no Indications of changes In the absorption or spectral Index during 
the more rapid changes (Mushotzky, Holt and Serlemltsos 1978; Baity et al . 
1983). 

There have been a few observations of rapid variability made with 
the Einstein Observtory. Tananbaum (1980) reported strong evidence for 
variability from NGC 6814 on a time scale of 6 hours with an indication of 
activity at shorter time scales. In the same article, Tananbaum reported a 
flux increase from 3C 273. The flux of 3C 273 Increased by 101 between two 
observations separated by ~ 12 hours. Tananbaum (1980) suggests that the 
change could have taken place In ~ 6000 sec on the basis of a single 
unpublished data point. 

Tananbaum et al. (1979) reported that the luminosity of the quasar 
OX 169 varied by 1.5 x 10^ erg/sec. The reported time scale of 6000 sec 
Is based on the difference between two data bins with rather large 
uncertainties. If one Includes this statistical error then one finds that 
the observed change In AL/At could be smaller by a factor of two. 

Matllsky, Schrader and Tananbaum (1982) have reported evidence for 
200 sec variability from the quasar 1525+227. Their published light curve 
consists of 8 data bins (200 sec each) of which 2 bins, near the end of the 
observation, are high. Also published were the arrival times of all 107 


photons collected. 

1.6 Variability At Other Wavelengths 
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The most spectacular observation of rapid variability has been 
reported by Wolstencrott et al. (1982). They observed the BL Lac OJ 287 at 
a wavelength of 1.25 pm. Their published light curve shows that the flux 
can vary by a factor 2.0 on a time scale as short as a minute. At a z of 
.306, and assuming H 0 * 100 km/sec/Mpc, the change In the luminosity (in 
the observed band only!) Is 10 44 ergs/sec. As Wolstencrott et al. point 
out, such a rapid change In AL Implies relativistic beaming. Beaming helps 
in several ways; photons (In the Earth's frame) have higher energies, 
timescales appear shorter, and the flux is not radiated uniformly 
over 4ir sterradlans. The combination of these effects could bring the 
value of a L/A t down to the theoretical maximum (Fabian and Rees 1979). 
Wolstencroft et al. saw rapid variability in only one observation (out of 
12) of OJ 287. In addition, they have observed other QSO's and BL lacs for 
many tens of hours, without seeing such changes. Thus rapid variability Is 
very rare. 

In the radio part of the spectrum, Kikuchl et al. (1973) reported a 
20% decrease In the 7.2 cm radio flux during an observation of OJ 287 which 
lasted 100 min. There was only a very small Increase In the flux of 
3C 147, which was observed before and after OJ 287 and served as a 
calibration source. Efanov et al. (1977) reported a 20% change over 4 
hours In the 22 GHz flux from 3C 273. The change over 24 hours amounted to 
33%. Coe et al. (1983) searched nine active galaxies for 10.7 GHz 
variability. They considered time scales ranging from 25 min to a few 
days. They reported no strong evidence for variability on a time scale of 
less than 1 day but did see weak evidence for day-to-day changes in the 
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flux from 3C 273 and QSO 0241+62 . 

Lawrence et al. (1981) observed NGC 4151 at optical wavelengths for 
several hours per night for a week. They reported that the nightly means 
were constant to 0.2 mag (2%) over the entire week. Lyutyi (1977) reported 

that NGC 4151 continuously shows flares with a rise time of 2-10 days and a 

decay time of 20-50 days. Thus Lawrence et al. seem to have observed NGC 

4151 during a quiet period. On shorter time scales Lawrence et al . 

reported a slnjle .15 mag (15$) dip lasting ~ 2 hours. Since the colors 
did not change during this time, they conclude that a high thin cloud could 
have caused the dip. With this exception they conclude that they did not 
observe any significant variability on time scales ranging from 10 sec to 
several days. 

1.7 Outline of Research Presented 

The main research presented in this thesis Is a very sensitive 
search for rapid X-ray variability. In order to make this search, one must 
be able to extract source variability from X-ray data. Chapter 2 describes 
a method that works when the data bins have unequal lengths. Given an 
observation of variability one then attempts to measure some sort of time 
scale. Methods to do so are presented In Chapter 3. These methods have 
been tested with numerical simulations. The final preliminary step Is to 
evaluate the performance of the detector. This Is done In Chapter 4 where 
the amount of residual noise from various sources Is estimated. 

Chapters 5, 6 and 7 form the bulk of the astrophyslcal results 
presented. In Chapter 5 It Is shown that the majority of active galactic 
nuclei studied do not vary on time scales of less than 1 day. The test was 
sensitive enough so that In an observation lasting < 12 hours, variability 
with a time scale of a day or more could be detected. The exception to the 
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above was NGC 6814 which showed factor of two changes on time scales as 
short as 100 sec, but with no detectable spectral changes. Chapter 7 
attempts to reconcile the rapid variability of NGC 6814 with the general 
lack of variability observed from other objects. 


11 


2. SEPARATING THE SIGNAL FROM THE NOISE 

2.1 Chapter Overview 

In this chapter I will describe a general method to separate true 
source variability from photon counting noise. The method works for the 
case where the data bins have unequal lengths. It Is also possible to 
derive an upper limit for variability If the count rate Is consistent with 
a constant. Next, I consider some simple models for variability and show 
the sensitivity of the method. It Is always possible to obtain a more 
sensitive upper limit If one has some a priori knowledge of the shape cf 
the light curve. Various distributions and their moments are described. 

In the final section the method was tested with numerical simulations. 

2.2 Comment on Weighted Quantities 

For X-ray astronomy a dominant source of noise Is due to Poisson 
counting statistics. If the true mean counting rate Is R and one 
Integrates over a time Interval of At, then counting statistics will 
introduce an uncertainty « = (RAt)^ 2 In the total number of photons 

r 

collected. Since the total number of photons C ■ Rat Increases linearly 
with time the percentage uncertainty In the rate decreases 
as At Increases. One does not know the trtoe rate R and so Is forced to 
estimate this number from the data. 

If the rate Is a slowly varying function of time then one can 
estimate the true rate at a given time by summing over a large number of 
bins N. Thus 

i N C. . 

<R> = hat f cl s mr 
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where Cf Is the observed number of counts in bin i and Ctot “ E i c i* 
total number of counts observed. Since At 1/N we find that the average 
Is Independent of the number of bins into which the data was divided. 

When one considers weighting the data more care Is needed. If, for 
every observation of Cj, there exists an Instrumental error then a 
maximum likelihood analysis (see Bevlngton 1979, p. 69) shows that the best 
estimate for the mean Is given by 

<R> W * ^WjCj/AtEjWj 

2 

where < > w will denote a weighted average and ■ 1/a^ . One Is tempted 
to base the weight for an Individual bin on the observed number of counts 
In that bln. In Appendix A, I show that for small constant mean rate this 
will generate an Incorrect result for the average (no matter how many bins 
are considered). This Is due to the fact tha* for small C.j, the square 
root of Cj Is a relatively poor estimator for the true uncertainty. 
Therefore when one weights data It is better to base the weights on the 
expected number of counts rather than the observed number. If is large 
then percentage difference between the expected and observed count rates Is 
small. However, In section 4.6, I give an example where the difference Is 
notlcable for = 100. 

2.3 Separating Source Variance From Photon Noise 

The simplest way to detect source variability is to 
calculate x assuming that the source flux is constant. The size 
of Itself provides a useful measure of the source variance (see, for 
example, Boldt et al. 1975). Also, ^ tables provide an estimate of how 
likely the variance In the observed signal Is due to chance alone. The 
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disadvantage of x 2 Is that It measures the amplitude of source variance 
relative to photon noise. Thus different experiments with different 

p 

effective areas will give different values of x for the same signal. In 
this section I show how one can separate the source variance from photon 
noise. 

If we assume that the source flux F Is constant but that the 
exposure Ef (effective area Integrated over time) can vary from bln to bln 
then x 2 is given by 

, (C, - © ( ] ! 

X - *1 - C — - 

E • 

B E.j (Fj " <F>)*" (2.3.1) 

where <C>^ » <f> Ej Is the expected number of counts for bln 1 and Cf Is 
the observed number. Since x 2 is effectively the weighted variance, It Is 
correct to divide by the expected numbe, of counts for that bln as was 
argued In the last section. Since Ej/<F> Is effectively the weight, we 
find that longer exposures get more weight. 

Since x 2 is related to the weighted variance with 

X 2 * <W> <o 2 > u( (2.3.2) 

V W 

where x 2 is reduced x 2 > <W> = z and W.j is the weight given to bln 
1. The total variance Is the sum of photon noise and any excess 
variance o 2 . If we solve for the excess variance we find 
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/<W> - < 0p 2 > w 


• <"p>w [4 /(<w> <£«> - n 

2 

If tiiere Is no excess variance, we know <x^> s 1. 
requl re 


(2.3,3) 

Thus for consistency we 


J P : 


<0 -> * <W> 


(2.3.4) 


Since neither <W> nor <o^> Is a function of we find that In general, 

P e 


x 2 - 1 

Ay 


<0 Pw s <W> 


(2.3.5) 


The final result will be the error In the background subtracted 
source flux (F s ). Let £ be the total exposure (effective area Integrated 
over time) then 

C - C R 

Fj. = — £-2. (2.3.6) 


where C Is the total count rate and C B the background rate. The 
statistical uncertainty In F s Is given by 

var(F 3 l (2.3.7) 

E 

where I have assumed that C B Is a known quantity. Now E = where 

n tot la the total area and e Is the efficiency (e “ 1 if source Is on 
axis), Thus 


A 
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, E F s + A tot ir Fg EF s tE Fg/e 
var(F s ) ■ — = — — « — 2 — 2 — - — 


F $ + Fg/e 


(2.3.8) 


where the background flux F B Is not a function of e. 

For the observations reported In this thesis the count rate Is 

dominated by a background which will be shown to be quite constant. For 

the smallest bln size used, 5,12 sec, the minimum number of counts per bln 

Is ~ 56. For this large number of counts we expect that the error made by 

setting a ^ (Instead of <C>) to be small. However, for the blank sky 

analysis (section 4.6) where I added together ~ 50000, 5,12 sec rates I 

find that using Introduces an unacceptable error If <C> < 120. If <C> 

exceeds 120 then It makes little difference which expression one uses 

for even for the large number of bins considered. 

2,4 Calculation Of Upper Limits 

If reduced x Is less than or close to 1.0 then statistical 

fluctuations will dominate any true source variability. For this case 

equation (2.3.5) breaks down, and It Is better to consider an upper 

limit. I will base my upper limits at 90% confidence level. Thus, a 90% 

, p 

limit means that, If there were no source variability, then x would exceed 
the given value 10% of the time due to chance alone. 

For a large number of degrees of freedom v we find that the 
probability of exceeding a given value of x 2 Is approximately P (X) where 



/2vT 


(Equation 26.4.11 of Abramowltz and Stegun, 1970) and P(X) Is given by 
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Table 26.1 of Abramowltz and Stegun. For large v, a value of X * 1.3 gives 
a 90% confidence level. What Is surprising Is that this value provides a 
good estimate of the 90% confidence level for all v. In fact X*1.3 
corresponds to 9.2% for v=l, 10.0 for v=2, 10.5% for v=10 and then drops 
back down to 10.0% by v=200. 

Since the equation for X Is so easy and accurate I have used It all 

the time, even in computer programs where one could easily consider higher 

order corrections. Therefore, a 90% confidence upper limit is obtained by 

? 2 

setting X * 1.3, solving for x and then using this value of x to obtain 
the excess variance. Thus the observed value of x 2 is not used to 
calculate upper limits. 

2.5 Sensitivity 

It is Important that one understand what the upper limits actually 
mean. Therefore, In this section I consider several possible models for 
variability and work out what minimum amplitude signal would, on the 
average, be detectable at the 90% confidence level. I point out that. If a 
source population did vary with this minimum signal, then the observed 
variability would exceed the upper limit 50% of the time and for the 
remaining observations only upper limits would be reported. This Is 
because, on the average, the observed excess variance will be above the 
true excess roughly half the time. 

Consider an observation which lasted for N bins. Further assume 
that N-n bins were at one level y L and the remaining n bins were higher by 
an additional amount 2a. Thus the variance due to this step function type 
change is given by 


°e "TTT + (N-n) (y L +2A-u) <J } 


(2.5.1) 
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It Is easy to show that p = y^+ 


Using this expression gives 



Now we use 


_ 4a 2 (N-n)n 

~ m: — r — 

4 = 1,3 4 and s ° ive f ° r 


(2.5.2) 


1/2 

24/c p - 1.34 (N-l) 1/4 


(2.5.3) 


where 2A/Op represents the total change measured In terms of the 
uncertainty for a single bln. 

For n = 1 we have the case of a single flare In the observation. 
For this case (2.5.3) reduces to 


2A /a ■ 1.36 N 1/4 (2.5.4) 

r 

which shows that as one considers more bins, the minimum detectable flare 
size, 2A/ap, Increases. This makes sense since the probability of seeing a 
given no fluctuation Is much higher as one considers more data. In Table 
2.1, I evaluated (2.5.4) for certain values of N. For 1000 bins a 
7.6 a fluctuation Is needed to give an unacceptable x for the entire data 

o 

stream. A 7.6 a bln contributes 58 to x and so our reasoning says that 
50£ of the time one can hide a Ax 2 of 58 In the remaining 999 bins. It Is 
unlikely that a 7.6 a fluctuation Is due to chance alone If one has only 
1000 bins. It Is important to realize that this method will not provide 
the smallest upper limits. If one Is Interested in single bln flares. 
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Table 2.1 

Magnitude of Change for 50% Detection 


One Bln High Half The Bins High 


N 

n 

2A/* p 

N 

n 

2A/o, 

10 

1 

2.4 

10 

5 

1.52 

100 

1 

4.3 

100 

50 

.86 

1000 

1 

7.6 

1000 

500 

.48 


Also In Table 2.1 I have listed the sensitivity for the case where 

half the data Is at p-A and the other half at p+A. We can see that as one 

considers more bins, one Is more sensitive to such a change. Although 

there Is no reason to do so one can assume that the first half Is at one 

value and the second half at the other. For the 1000 bln case, a 

.47 cr fluctuation would be detected ~ 50% of the time. However a step 
P 

o 

function model would be expected to drop x by 58 which would be very 
significant. In general, If one considers a specific model, then one can 
obtain limits more sensitive than the value obtained In Section 2.4. 

2.6 Autoregressive Process 

The autoregressive process Is a simple method used to generate data 
with properties similar to observed data streams. Since It Is simple, one 
can quickly determine the parameters and then proceed to generate a lot of 
phantom data which looks similar to observed data. This allows one to 
estimate the uncertainties In the derived parameters. 

The n^ order autoregressive process (see Jenkins and Watts 1968) Is 


I 
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defined for a discrete set of points by 


Y 


1 


N 

* E co 

J-l J 




(2.6.1) 


where are the parameters and 7, is a random variable. Scargle (1981) has 
written ari excellent introduction to time series analysis. He shows the 
relationships between an autoregressive process and others, such as a 
moving average modal and a shot noise model. 

Consider the first order process 


Y i ■ °Vi + z i 


To determine the average value for Y consider 


( 2 . 6 . 2 ) 


<Y.j> = a <Y i>4 > + <Zj> 


(2.6.3) 


Now if we assume the process Is stationary, then <Y^> = <Y-j_j> s <Y>. If 
we define <Z> = <Z-j> : we then find 


<Y> = St. 
1 T=a 


(2.6.4) 


To determine the variance of Y one has to consider 


<Y 1 Y f > = <(c*Y i „ 1 + Z 1 )(aY 1 _ 1 +Z 1 )> 


<Y 2 > = a 2 <Y 2 > •+ 2a <YZ> + <Z 2 > 


(2.6.5) 


ft 
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Since Zj and Y^i are uncorrelated* <YZ> ■ <Y><Z>. It Is possible to 
substitute (2.6.4) Into (2.6.5). When one solves for the variance 
<(Y-<Y>) 2 > one obtains 


var(Y) » . var . ( , z ) (2.6.6) 

1 - c* 

Finally, It Is possible to show that the skewness or <(Y-<Y>) 3 > Is given by 


skew(Y) = — (2.6.7) 

1 - a 

Equations (2.6.6) and (2.6.7) Indicate that the shape of the distribution 
of Y Is given by the shape of the Z distribution but the scaling from Z to 
Y Is controlled by the value of a. 

2.7 Moments 

In this section the moments of various distributions are tabulated 
for easy reference. 

First consider a distribution in which the numbers are confined to a 
range utA. A one parameter model for such a distribution could be 


P n (x) = tfk {1 * 


where |x| < A. Notice that as n Increases the 11m P n (x) = which 
corresponds to a uniform flat distribution. The variance of P n Is given by 


var(P n ) = / P n (x) x^dx 
-A 


a 2 n+1 

riT 
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As n ranges from 1 (triangular distribution) to » (uniform distribution) , 
the variance changes by a factor of 2. Although the variance Is a 
continuous function of n, the dependence on n Is quite weak. Finally the 
skewness of P n Is always zero since skewness measures the asymmetry of the 
distribution. 

A Gaussian distribution is given by 


P G (x) 


. _L_ 

a/Z n 


e 



p 

The first three moments are y,a ,0. Again we have a symmetric distribution 
and hence the skewness Is zero. 

The Poisson distribution, 


PpfK.nl 




I 


Is Interesting In that the discrete case Is easier to manipulate than the 
continuous case. For exammple consider computing the mean 


<x> = z 
x=o 


* xT e " M = 


Tx^TT! 




It Is now possible to set y=x-l and to do the summation to obtain <x> = y. 
In general It Is possible to express the n tfl moment In terms of (n-l)^ h 
moments. Thus one can slowly proceed to higher moments finding 

t 

that var(P ) = y and skew(P ) = y. Before the reader jumps to a hasty 

r r 

generalization let me point out that the 4 t * 1 moment Is 3y 2 +y. 
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Astrophysicists often use distributions which are generated by shot 
noise models. In such models one assumes that the observed light curve Is 
composed of a large number of superimposed flares. In the standard shot 
noise model one assumes that all the flares are identical. Let H(x) be the 
intensity of a single flare after time t from the beginning. Then the 
observed flux Is 

I(x Q ) c / Z(t) H(x q -x) dx 

where Z(t) = 1 when a flare begins and 0 otherwise. If ** denotes 
convolution then the above can be written as 


I * Z**H. 


where one thinks of H as being the Impulse response function and Z as a 
series of Impulses. 

To compute the mean we assume that the process Is stationary (the 
mean Is not a function of time). Since It makes no difference when the 
mean is computed we will do It at time t = 0. Thus 

<I> = <I(0)> = </ Z(x.) H ( -t ) dx> 

Since the expectation Is calculated for a given time we can Interchange the 
order of integration and expectation to obtain 


<I> = / <Z(x) H ( —x ) > dx 
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* / <Z(t)> H(-t) dt 

since H ( t) Is a constant for any given t. Since the process Is stationary 
the <Z(t}> Is constant and Is generally associated with the symbol X. We 
note that x Is the mean shot rate and has units of Inverse time. So at 
last 

<I> ■ X / H{ t) dx. 

The variance Is slightly more tricky. We have 
var(I) * var( 1(0)) 

a < / Z(x) H(x) dx / Z(x') H(x') dx'> 

= //< Z{ x ) Z(x'}> H(x) H(x') dxdx' jj 

j , 

Now we must make use of the assumption that the Z's are uncorrelated. jt 

II 

Therefore j 

<Z{ x ) ZX (x')> = 0 If x*x' 

If x s x‘ we have <Z z (x)> but, since by assumption Z Is either 0 or 1 we .i 

9 \ ' 

have <Z (x)> = <Z(x)> = X. Thus m 


var(I) = X / H 2 (x) dx. 



The same reasoning can be applied to the skewness to obtain 


skew(I) « x J H 3 (t) dx 


Table 2.2 


Moments of Various Distributions 


DISTRIBUTION 

MEAN 

VAR 

SKEW 

n+1 M ixl n v 
TnE (1 " nA } 

P 

A 2 {n+1) 

TTnTTT 

0 

Gaussian 

p 

2 

0 

0 

Poisson 

P 

p 

P 

shot noise 

Xxh 

u2 

Xx-jj- 

h 3 

Xxj 

H(t) = he‘ t/x 




shot noise 

Xhx 

Xxh 2 

Xxh 3 

H(t) = h, t < x 




AR 

<z> 

var(Z) 

skew(Z) 


1-a 2 

“ — r 

1 - a 


AR (uniform Z) 

<Y> 

<Y> 2 1-ct 
“3 TPa 

0 
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2,8 Numerical Results 

Analytical results should be complemented with numerical work, since 
using both methods often allows one to reach greater understanding. In 
this section I briefly discuss ny program {called PHANTOM) which generates 
"phantom" data. This program has been used to test other programs as well 
as simple (analytic) models. 

PHANTOM produces a data file In exactly the same format as .typical 
X-ray rates files. Thus programs can use phantom data without 
modification. It Is currently possible to generate two types of time 
series: an autoregressive process (first or second order) and shot noise 
(either square or exponential shots). To the basic time series one can add 
a constant term, sine waves, or noise. The program first generates a 
series of "true" rates p^. If the user requests It the program then 
constructs a series Xj such that x-j has a Poisson distribution about p^. 

The method used allows the user to generate the same series of (for the 
same starting seed) whether or not noise is added. The noise term Is 
Poisson distributed If there are fewer than 50 counts per bln. For larger 
rates the FORTRAN code which generates Poisson noise falls and so a 
Gaussian approximation Is used with a 2 a p. Although the difference 
between Poisson and Gaussian noise Is small for large rates, I find that 
the difference Is easily detected (see Table 2.3). 

As a demonstration of what one can discover consider Table 2.3. On 
the first line, I consider the mean, variance and skewness of a pure 
autoregressive (AR) process. The expected values are 100.0, 175;4, and 0.0 
respectively. The plus or minus values are the single sample uncertainty, 
as determined by running the program 10 times. Notice that as one goes to 
higher moments the uncertainty Increases drastically. 
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On the second line of Table 2.3 I give the observed mean, variance 
and skewness for constant rate but with added noise. The expected values 
are 100, 100 and 0.0. If the noise had a Poisson distribution then the 
skewness would be 100 but as I mentioned I am using a Gaussian 
approximation and so skewness * 0.0. Again the uncertainty Increases with 
the higher moments. 

Finally, In the third line of Table 2.3 I consider adding noise to 
the autoregressive process generated In the first line. Notice that the 
variance of the sum Is equal to the sum of the variances due to the AR 
process and due to the noise. However the skewness does not add. To see 
this let P(p) be the probability that the signal Is at a level 
of p (without noise). Further let P(x; p) be the probability of seeing x 
photons If the expected number Is p. The probability of seeing x photons 
during the entire observation Is 

P(x) s / P(x; p) P(p) dp 

The expectation of x n Is given by 

<x n > = // x n P(x; p) dx P(p) dp. 

If It Is possible to do the Inner Integral, then one can express <x n > In 
terms of the various moments of p(p). For Gaussian noise we find that 
skew(AR+no1se) * skew(AR) + 3var(AR). For Poisson noise (see Appendix B of 
Sutherland, Welsskopf and Kahn 1978), we find skew(AR+no1se) = skew(AR) + 
3var(AR) + mean(AR). 
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TABLE 2.3 

MOMENTS OF VARIOUS DISTRIBUTIONS {EXPERIMENTAL) 



MEAN 

VAR 

SKEW 

pure AR 

100.0*1 .3 

171 .7±16 .3 

71 • 1±307 .8 

pure noise 

99. 6± .2 

100. 9* 2.1 

-5.3* 72.3 

AR+nolse 

99.6*1 .4 

275.1*18.7 

594.1*444.8 
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3. DETERMINATION OF TIME SCALES 

3.1 Chapter Overview 

O 

In the last chapter we saw how to use x first to detect variability 
and second to determine the source variance. Given an excess variance one 
next proceeds to measure the time scale for variability. There Is no 
standard method that works for all cases, so some human judgement Is often 
used. If the observation lerjth Is short, relative to any apparent time 
scale, then one can measure the time for the source to double Its flux, the 
source doubling time. For a longer observation the peak-to-peak time, or 
the time Interval between peaks In the flux, Is perhaps the physically most 
useful time scale. For a long observation one can apply powerful 
statistical methods to extract the time scale. These methods Include 
computing the auto-correlation function (related to the power spectrum) and 
computing variance as a function of the bln size. 

3.2 Short Observations 

If the apparent time scale is much longer than the observation 
length then one will only see slow trends In the data. In order to 
characterize the data one generally fits a polynomial to the observed 
rate. If a straight line provides an acceptable fit to the data then the 
"doubling time" Is a useful quantity. I define the doubling time to be 
the mean Intensity divided by the slope. Thus the doubling time measures 
how long It, would take the source to double Its flux if the observed trends 
were to continue. 

The doubling time Is useful for several reasons. First, flickering 
sources rarely vary by as much as a factor of 10, and so the time for the 
flux to double will provide an order of magnitude estimate of the true time 
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scale. Secondly, the doubling time prices an estimate for how long to 
observe the source when requesting additional observations. 

If a slow trend fits the observation then the slope of the trend Is 
Interesting since Fabian and Rees (1979) (see Section 1.4) have proposed a 
theoretical upper limit to AL/At. For variability near this limit, one can 
strongly constrain various physical models (see Chapter 6). 

3.3 Medium Length 0. xervatlons 

When one nas observed a source over only a few variations It Is 
difficult to define the observed time scale in a formal manner. In this 
case one can examine the light curve by eye. The time scale that one looks 
for is the time from one extreme to another for Individual events. Due to 
noise It Is not generally possible to locate extrema precisely. However 
this peak-to-peak time scale Is physically very Important since It provides 
an upper limit to the source size. Thus in the absence of relativistic 
bulk motion 


At < R/c 
PP ~ 


This Is due to the fact that some parts of the source are closer to the 
observer than others, and can be separated by at most a distance R along 
the observer's line of sight. 

3.4 Long Observations: The Auto-correlation Function 

When the observation lengtn is long relative to the time scale of 
variability then one n compute the auto-correlation function. Given a 
set of observations with statistical errors we compute the 
auto-correlation function with 
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p(u) 

r(u) 


= r(u)/r(0) 

1 N'“ < Y 1 - <Y> w )(Y t+u - <Y V 


(3.4.1) 

(3.4.2) 


The variable u Is the lag time. Jenkins and Watts (1968 p. 772f) showed 
that using 1/N Instead of 1/ ( N-u) In the definition of r(u) Introduces a 
bias, however the biased function has a smaller mean square error. 

The definition of p(u) shows that It Is directly related to 
reduced x 2 . Thus In order to compute a meaningful auto-correlation 
function, x 2 should be significantly different than 1.0. Some people feel 
that the auto-correlation function can magically pull a signal out of quiet 
data. Although this may be possible, It has never happened In any of my 
tests. The auto-correlation function Is useful only If a x 2 test shows 
that a constant Is an unacceptable fit to the data. 

Welsskopf, Kahn and Sutherland (1975) pointed out that photon noise, 
being uncorrelated, contributes to p ( 0 ) but not to p(u) If u Is non-zero. 
They recommend correcting for noise by multiplying p(u) by a correction 
factor F for all non-zero lags. An expression for F Is 


F *» 



(3.4.3) 


where N data points have been used. This expression makes It clear that 
If x 2 ~ 1.0 then the correction factor will be quite large, l.e. the ratio 
of sign?) to noise will be very small. 

The auto-correlation function of shot noise Is given by 

I 

I 


i \ 

s 


1 



t 


r(u) = \ / H(t) H(t+u) dt 


(3.4.4) 
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We see that the auto-correlation function Is no more than the convolution 
of the Impulse response function with Itself. If H(t) * he" t//r then 

r(u) * e" u ^ T (3.4.5) 

and of course 

p(u) * e“ u/T (3.4,6) 

For a first order autoregressive process the autocorrelation function Is 

p(u) = a u = e" u/T (3.4.7) 

where t ** -1.0/ln(a) (see Jenkins and Watts 1968). This clearly shows the 
relationship between a and x In addition to the similarity between shot 
noise and an autoregressive process. 

The maximum likelihood estimate for a Is p(l) (see Jenkins and Watts 
1968). Although one feels that it should be possible to obtain a better 
result using more than one point, numerical tests have shown that the best 
estimate is p(l) alone. This Is due to the correlations in the data, which 
cause the deviation of p{2) to be In the same direction and slightly 
greater than the deviation of p(l). For no Poisson noise it is possible to 
estimate the uncertainty In p (1 ) with 

< 7 2 . [1 - p(l) 2 ]/N 


(3.4.8) 
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(see Box and Jenkins 1976, p. 34). 

The above were demonstrated numerically. A first order 
autoregressive process was used to generate 2000 points. The random 
variable Z varied uniformly between 0 and 20, Since a * .9 the mean rate 
was 100 ct/bln. Ten runs were made without adding counting noise and the 
mean for these runs was .8970 (see Table 3,1). The deviations about the 
mean were ± .0093, compared to a predicted quantity of .0097 using equation 
(3.4.8). When photon noise was added the mean value for p(l) dropped to 
.5731. However when one corrects each p (1 ) for noise, the mean rises to 
,8936 which again Is very close to .9, the expected value. 

3.5 Long Observations: Yarlance vs. Bln Size 

Another method of estimating the time scale for a long observation 
Is to compute the excess variance (section 2,3) as a function of bln 
size. This method works because, for long bln lengths (sizes), one Is 
averaging over several fluctuations. This causes the variance to 
decrease. As the bln size decreases then, at some point, all fine scale 
structure Is "resolved out". Thus for short bln sizes the variance is 
constant. The overall effect Is that, as one Increases the bln length, the 
variance Is constant for a while, begins to roll over when the bln length 
equals the characteristic time scale for variability and for large bln 
lengths decreases monotontcally. 

Define Bn to be the binned data where n data points are averaged 
Into a single bln. Thus 


i n- 
" n * =0 V (+j 


(3.5.1) 
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TABLE 3.1 

AUTOCORRELATION FUNCTION OF FIRST ORDER AUTOREGRESSIVE PROCESS 



No Photon Noise 

With 

Photon Noise 

Trial 

p (1) 

p(l) 

F *p ( 1 ) 

1 

.8868 

.5229 

.8644 

2 

.9129 

.6495 

.9535 

3 

.8922 

.5698 

.8980 

4 

.9008 

.5782 

.8835 

5 

.9007 

.5588 

.8924 

6 

.8966 

.5947 

.9301 

7 

.9053 

.5654 

.8413 

8 

.8801 

.5339 

.8751 

9 

.8951 

.5711 

.8863 

10 

.8991 

.5871 

.9118 

x » 

.8970 

.5731 

.8936 

= 

.0093 

.0374 

.0322 
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In the remainder of this section I will suppress the subscript 1. For a 
stationary process we find 

i n-1 

<Bn> ■ ■=■ Z <Yj> * <Y> (3.5.2) 

n J«0 J 

Now consider the variance of the binned signal, 

var(Bn) = <{E * Y.) (S 1 Y. )> - <Y> 2 (3.5.3) 

rr j*0 J k=0 K 

To evaluate the product of the sums, It Is Important to realize that the 
expectation of the product of any two terms Is a function of only the 
difference between the Indices. Thus <Y^ + ^ Y^> = <Y-|+2 Yi+l* etc * In 
product of the two sums there will be n terms with the same index. If the 
Indices are not the same then there will be 2(n-j) terms In which the 
indices differ by j. Plugging this Information Into (3.5.3) gives 

n 1 

var(Bn) = 4 <Y 2 > - <Y> 2 + E (n-1) <Y Y.> (3.5.4) 

n n i 

where <Y Y^> Is the expectation of the product of two terms with Indices 
differing by an amount 1. From the defl nation of the auto-covariance 
function we have 

r(1 ) = <(Y - <Y>) (Y i - <Y> )> 

= <Y Y^> - <y> 2 (3.5.5) 


Substituting (3.5.5) Into (3.5.4) gives 
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n-1 


var(Bn) « I <y 2 > - <Y> 2 + is (n-1)(r(1 )+<Y> 2 ) 


n 1=1 


(3.5.6) 


Now 


Kl 1 (n-1) <Y> 2 = (1 +i) <Y> 2 
n 2 1=1 n 


so (3.5.6) reduces to 


(3.5.7) 


var(Bn)/var(Y) = I + L. S’ 1 (n-i) p(1) (3.5.8) 

n n^ 1=1 

where I have used the definition of the auto-correlation function, p(1). 

3.6 Examples of Variance vs. Bln Size 

First assume that there are no correlations In the data (white 
noise). For this case (3.5.8) reduces to 

var(Bn)/var(Y) = l./n (3.6.1) 

which shows that the standard deviation Is proportional to l.//rf which Is 
what one expects for Independent samples of a distribution. In Figure 3.1 
the curve a = .0 corresponds to the case of no correlations In the data. 

For a first order autoregressive process (p(1) = a*). It Is possible 
to do the summation to obtain 

var(Bn)/var(Y) = Lv - “ r - 2 ^ 1 -- QC -i (3.6.2) 

" TT ^ r n 2 (l-«) 2 
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In Figure 3.1, I have plotted the square root of (3.6.2) vs. log 2 (n) 
for a = .50, .90 and .99. Since the time scale Is ~ -1/1 n{ a) we see that 
the curve remains horizontal out to large blnnlngs If the time scale Is 
rather long. 

Finally the program to separate the excess variance from a data 
stream was tested. The results are shown In Figure 3.2. Since the data 
stream consisted of only 2000 points, the variance for the largest bln size 
was computed with only two binned points. Therefore the agreement with the 
theory Is very good. 
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FIGURE CAPTIONS 

Figure 3.1 — The square root of (3.5.8) plotted as a function of n for a 
first order autoregressive process. For a * 0, no correlations 
exist In the data whereas for a ** .99 the characteristic time scale 
is 99.5 bins {log 2 (99.5) - 6.64). 

Figure 3.2 — The results of a numerical experiment. Ten data streams were 
generated with a first order autoregressive model. For each data 
stream a = .90 and 2000 points were generated. The curve represents 
the average and the error bars represent the single sample 
uncertainty. 
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4. IN FLIGHT PERFORMANCE 

4.1 Chapter Overview 

In this chapter I will discuss the in-flight performance of the 
Goddard instrument on board HEAO-1 and how this performance affects timing 
experiments. The Goddard experiment rejects particle events with over 99% 
efficiency; this small residual background noise can be further reduced by 
selecting times corresponding to low values of Mcllwain L. Detector 
offsets were refined in flight and any remaining offset uncertainties can 
be neglected. The uncertainties in the aspect solution are also 
discussed. Finally, a variance analysis was performed on blank sky 
points. These observations prove that the systematlcs are quite small. 

4.2 The Instrument 

The Goddard experiment on HEAO-1 was a gas f 111ed proportional 
counter. For a complete discussion of the detector see Rothschild et al. 
1979. Here I will only be concerned with some of the general properties of 
the Instrument. X-rays enter the gas volume via collimators on the front 
of the Instrument (see Figure 4.1). These X-rays photolonlze the gas 
(either argon for the medium energy detecter MED or xenon for the high 
energy detectors HED's). These photoelectrons produce a local region of 
ionization which gives rise to a current pulse. In Figure 4.2 a cross 
section of the instrument is displayed with the anode and cathode wires 
coming out of the page. The voltage Is controlled such that the current 
generated Is proportional to the Initial energy of the X-ray. 

Proportional counters are relatively Inexpensive and can be built 
with large collecting areas. In Appendix E, I point out the importance of 
large area detectors. Proportional counters (like most X-ray detectors) 
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also respond to the passage of charged particles. In the near Earth 
environment the count rate due to charge particles can be 100 times greater 
than the X-ray count rate for extragalactlc sources. This background Is 
also a strong function of the spacecraft's location In the Earth's magnetic 
field. Thus, as the spacecraft moves along Its orbit, one expects a time 
variable background signal. 

Charge particles rejection makes use of the fact that X-ray- events 
tend to be localized in the volume whereas particles tend to Ionize along 
tracks (see Mason and Culhane 1983). Methods of rejection fall Into two 
broad classes; pulse shape discrimination and anticoincidence. It should 
be clear that localized X-ray events should give rise to pulses with a mean 
shape different than the pulses produced by particle tracks. This method 
has been successfully used by others (see Gorenstein and M'WiewIcz 1968). 

The Goddard Instrument did not use pulse shape discrimination but 
rather relied entirely on anticoincidence. To use an anticoincidence 
method the gas volume Is divided into many small cells (see Figure 4.2). 

As Implied by the figure the cells are long and thin. Thus It Is possible 
for a particle travelling parallel to a wire to stay entirely within one 
cell and not trigger the anticoincidence logic. Due to the small size of 
the cells this should be quite rare (~ .07%). Newer detectors (Bally, 

Smith and Turner 1978) have an "end veto" system to detect particles 
travelling along a wire. 

Low energy particles can also escape the anticoincidence logic, by 
losing most of their energy to a single cell. To provide protection from 
these events an additional surface layer is provided which Is not exposed 
to X-rays. For MED and HED2 three sides of the gas volume are protected by 
such a system (the V2 rate counts triggers on the sides and VI on the 
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bottom). HED1 and HED3 are protected along the top surface by a veto 
layer. In order not to attenuate the X-rays, the gas In the veto layer Is 
made of a low Z material (propane). The propane does respond to parHcles 
and so provides the veto signal. 

Finally It Is possible to monitor the Internal background of the 
Goddard detectors. Each detector has two different size collimators, one 
having roughly twice the field of view of the other. Thus any Isotropic 
emission, such as X-rays from the X-ray background, that comes down the 
collimator Is detected twice as strongly In the large field of view. The 
count rate from non-coil Imated events will be the same In both fields. 

4.3 Particle Discrimination 

For the Goddard Instrument the total number of anti coincident events 
are recorded every 40.96 sec. Each 40.56 sec Interval Is called a major 
frame. The anticoincident rate Is proportional to the number of charged 
particles transverslng the detector. To test for charged particle 
contamination I have plotted the total xenon detector (HED3) count rate vs. 
anticoincidence In Figure 4.3. The data were taken from a blank sky 
observation (see section 4.6) on day 632 of 1977 (Sept. 24, 1978). Most of 
the scatter in the plot Is due to counting statistics. However, a slight 
upward trend can be seen. A least squares straight line was fitted to the 
data and the fit parameters are listed In Table 4.1 The slope of the line 
Illustrated In Figure 4.3 indicates that the variable rate Is at a level of 
only .59% of the anticoincidence (particle) rate. Also listed In Table 4.1 
are the fits for various discovery scalers (DS). Discovery scalers 1 and 2 
are the rates for the two fields of view for the first layer only of the 
detector. Discovery scalers 3 and 4 are the second layer rates. Discovery 
scales 1 and 3 always corresponds to the rate In the 3° field of view. 
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TABLE 4,1 

SLOPE OF BEST FIT LINE 




kate vs. 

Rate vs. 

Dlt 

DS 

Anti (percent) 

Mcllwaln L ct/(sec-L) 

HED2 

Total 

1.46 ± .18 

4.5 t .8 


1 

.62 ± .10 

2.4 i .4 


2 

.59 ± .12 

1.8 ± .5 


3 

.21 ± .04 

.3 ± .2 


4 

.04 ± .05 

.0 ± .2 

MED 

Total 

.78 ± .08 

2.8 ± .4 


1 

.36 ± .05 

1.3 ± .2 


2 

.35 ± .04 

1.1 ± .2 


3 

•04 i .03 

.1 ± .2 


4 

.03 ± .02 

.2 ± .1 

HED3 

Total 

.59 ± .04 

2.8 ± .2 


1 

.16 ± .03 

.8 ± .1 


2 

.21 ± .02 

.9 ± .1 


3 

.12 ± .01 

.5 ± .1 


4 

.10 i .01 

.5 ± .1 


! 

I 
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Table 4.2 


X 2 /DOF FOR VARIOUS MODELS 


DET 

DS 

CONSTANT 

ANTI 

MCILWAIN L 

CONSTANT 
(L < 1.2) 

HED2 

Total 

381.6/261 

305.4/260 

337.7/260 

234.6/197 


1 

320.5/261 

275.9/260 

281.6/260 

195.4/197 


2 

263.6/261 

238.8/260 

249.8/260 

187.5/197 


3 

260.2/261 

239.1/260 

257.1/260 

179.5/197 


4 

270.1/261 

269.6/260 

270.1/260 

200.3/197 

MED 

Total 

433.8/313 

327.4/312 

367.8/312 

243.5/233 


1 

378.7/313 

327.1/312 

344.9/312 

245.1/233 


2 

385.9/313 

300.2/312 

342.1/312 

239.5/233 


3 

292.7/313 

291.4/312 

292.0/312 

227.6/233 


4 

281.1/313 

280.1/312 

277.9/312 

210.9/233 

HED3 

Total 

822.3/484 

533.9/483 

597.7/483 

368.4/299 


1 

572.6/484 

528.5/483 

530.9/483 

339.4/299 


2 

620.5/484 

497.1/483 

534.3/483 

320-5/299 


3 

599.7/484 

511.1/483 

543.5/483 

365.1/299 


4 

645.7/484 

560.5/483 

572.6/483 

370.2/299 
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Several other points can be seen from Table 4.1. First the 
contamination In both fields of view Is always about the same, which 
Indicates that the particles not being vetoed are not coming down the 
collimators. This also shows that when one considers the "true" X-ray rate 
by subtracting the small field of view from the large, the residual 
contamination will effectively be zero. Another effect Is the 
contamination In the first layer of HED3 Is only one-third the level for 
HF.D2. Except for the propane veto layer on HED3 both detectors are very 
similar. This shows that the propane layer reduces the background noise by 
a factor of 3. Finally, the second layer of the MED (which Is twice the 
thickness of the second layer of HED3) has an extremely small 
contamination. 

In Figure 4.4 I have plotted the anticoincident rate vs. the 

Mcllwaln L value (see Appendix B). It Is clear that the two rates are 

correlated. In Figure 4.5 I have plotted HED3 total rate vs. L. Again a 

slight upward trend can be seen In the data. The fit parameters are listed 

In the last column of Table 4.1. A good approximation to the slope for 

2 

half the detector Is 1 ct/sec per unit change of L. Comparing x (Table 
4.2) for the fits to Mcllwaln L and the anticoincidence rate shows that 
anticoincidence is a slightly better model. However, we are not concerned 
about an exact fit at this point. 

It appears from Figure 4.5 that the experiment spends much time at 
low values of L with low anticoincidence rates. In Figure 4.6 I display a 
histogram of how often various values of Mcllwaln L occur for twd separate 
pointed observations. These histograms clearly show that excursions Into 
large values of L are quite rare. It therefore was decided to discard data 
at large values of Mcllwaln L. Thus for the remaining of this work I only 
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consider data with Mcllwain L values of less than 1.2. This means 
that ~ 1/3 of the previously good data are discarded. In Table 4.2 I 
show x 2 for a constant for the data with L < 1.2. As can be seen the 
amount of data was reduced by a factor of 1/3 and the excess variance was 
reduced by a factor of three. 

4.4 Detector Offsets 

The total number of source counts observed is related both to source 
flux F (counts per area per time) and also to the effective area A. Thus 

C s = FA At 


For the Goddard detector we have 


» - 11 - 


Where A^ ~ 400 cm, 8 Q is distance in degrees from the center of the field 
of view to the source measured in the off-scan direction. The 
quantity e g is the distance measured in the along scan direction. 

Thus (0 a * 0 o ) form a coordinate pair in an orthonormal system. The 
angle e m is the opening angle in the along scan direction and can be 1.5, 
3.0° or 6.0° depending on the detector. The opening angle in the off scan 
direction was always 3°. The above shows that the detector has a simple 
pyramid response to off axis sources. 

It would be nice if one could measure (e, ,e ) directly. However 

d 0 

what is measured is the look position of the star trackers. Therefore one 
must determine a set of rotation angles (small) which will give the look 
position of the detector and this must be done in flight. Jean Swank and I 
determined these angles by looking at the area corrected rates as the 
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spacecraft scanned over the Crab nebula. If the rates are 'U' shaped the 
opening angle needs to be changed; If the curves are "S" shaped, the 
offsets need changing. In practice this procedure requires that one 
guesses new parameters and iterate a few times. In Figure 4.7 I have 
Illustrated the rate for one scan over the Crab showing a flat response. 

4.5 Pointed Aspect Errors 

For the pointed data used In this thesis an additional aspect error 
existed. Th* was caused by the fact that HEAO-1 was designed to be a 
scanning spacecraft. Therefore, the standard (NASA) software package to 
determine aspect information from star tracker data, only worked while the 
spacecraft scanned. During pointed observations the aspect was 
extrapolated from the last scanning fix using only gyroscope data. The A-3 
experimental group did work out Improved solutions based on the star 
tracker data which they were to provide to us. 

We obtained ~ 25 aspect solutions from the A-3 experimenters. We 
found that for most cases the Goddard-supplied aspect was accurate to 
better than .05°. An error of this size would produce an additional noise 
In ctj/I of 2 % for our 3° x 3° field of view. This signal can only be 
detected for the brighter soures. Since we used the A-3 solution for most 
of the brighter sources and for the more variable weaker ones, we believe 
that the residual variations found for some of these sources are not 
related to aspect uncertainties. 

4.6 81ank Sky Analysis 

HEAO-1 had instruments looking out in two opposite directions which 
is acceptable as long as the spacecraft scans. In general, during a 
pointed observation, only one side could be targeted for a source. Thus 
during a mlnus-y point, when the other side of the spacecraft was pointed 
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at a source, the Goddard Instrument was effectively pointed at a random 
place on the sky. The rates for all mlnus-y points were examined. If the 
field contained a known source then the rates were only quickly examined. 
This study turned up a solar flare that produced a flux In our detectors 
(see Appendix C). Also discovered was a flare that appears to be from a 
new "nearby” galactic source (mentioned In Appendix C and discussed In more 
detail In Tennant and Swank 1983). 

If there were no HEAO-1 sources (flux > 2 x 10"^ erg/cm^-sec) In 
the field of view then the point was flagged as a "blank sky" 
observation. There were 49 points which were selected for the blank sky 
analysis. Basically one determines the excess variance for each 
observation in a manner similar to the active galaxy study. This analysis 
will provide an estimate for the residual non-source noise. In the next 
chapter I present a histogram of Og/B for the points separately, where one 
can compare this histogram to a similar one made for the active galaxies. 

Here I will consider the results for all the blank sky observations 

2 

added together. To do this one adds together the observed x values for 
each point separately. Since the background Is a function of location on 
the sky, a separate background was determined for each pointed 
observation. The background rate was assumed to be constant throughout the 
point. 

In Figure 4.8 I display o B /B (in percent) vs. bin size for bin 
sizes ranging from 5 sec to 90 min. As can be seen from the figure there 
is excess noise at the 1.2% level for the xenon detector and 1.6% level for 
the argon detector. Since the background counting rate In the xenon 
detector Is ~ 1.4 times larger than the corresponding rate for the argon 
detector, we find the same, absolute amount of excess noise in both 
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detectors. Figure 4.8 also shows that the excess variance drops for the 
largest bin sizes of roughly 90 min. This indicates that the noise has a 
time scale of ~ 90 min or one spacecraft orbit. 

In searching for the source of the excess variance one has two 
clues. First, the signal has the same magnitude in both the xenon and 
argon detectors and second the time scale is roughly 90 min. Although one 
might think that 90 min indicates residual particle contamination, I point 
out that this is not the case. All the 90 min time scale indicates is an 
orbital dependent effect, of which particle contamination is one 
possibility. I will argue that such contamination is not a problem. 

In section 4.3, I shewed that the count rate was proportional to 
Mcllwain L. Since we have selected only L < 1.2 we can estimate the 
residual contamination. If we assume a uniform distribution then from 
section 2.6 we find 

a 2 = j- = y (2 x .1 ct/sec) 2 

a = .06 ct/sec 

The factor of 2 appears since I am considering the entire detector and L 
ranging over ± .1 causes a flux change of t .1 ct/sec in half the 
detector. Since the Mcllwain distribution is not uniform but peaked we 
expect the actual contamination to be less than this estimate. The 
residual noise from the blank sky points has a level of 1.2% of 14 ct/sec 
or .17 ct/sec. Thus the Mcllwain L effect can contribute at most a third 
of the observed excess noise. 

Under some conditions we might expect the residual noise to be 
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proportional to the Internal background. In Table 4.3 I have broken the 
total detector count rate up Into an internal background and an external 
background (mostly X-rays from the sky). As can be seen the internal 
background for the xenon detector Is three times greater than the internal 
background for the argon detectors. However the external background 
contributes roughly the same rate in both detectors. Since the residual 
noise Is the same for both detector;, fluctuations in the Internal 
background could not give rise to the observed effect. 

To understand how the external background could give rise to such a 
signal one needs to know two things. First, as the detector is pointed at 
different areas of the sky the background varies by about 3%. Most of this 
variation is due to weak unresolved sources. Second, during a pointed 
observation the detector look position was not held fixed but was allowed 
to vary within a "dead band" which extends up to 1/2 degree from the target 
position. For typical observations the spacecraft tended to spend much 
time near the edge of the dead band and the look position tends to circle 
the dead band about once every orbit or 90 min. These two points lead me 
to believe that it is the spatial fluctuations in the X-ray background 
which give rise to the observed ~ 1 % signal with roughly a 90 min time 


scale. 
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TABLE 4.3 

TOTAL DETECTOR BACKGROUND COUNTING RATE 


Argon 

Internal 

1.9 ct/sec 


External 

9.4 ct/sec 


Xenon 

Internal 

5.5 ct/sec 


External 

8.7 ct/sec 
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FIGURE CAPTIONS 

Figure 4.1 — Cross sectional view of the xenon detector HED3 {or HED1 ) . 

Figure 4.2 (top) — Detector grid connections for the argon detector MED 
(also LED and HED2). X-rays from the 3° x 3° field of view are 
counted In the LI or L2 rates. Rates R1 and R2 will contain X-rays 
from the other field of view (3° x 1.5° for MED). Veto rates VI and 
Y2 shield the detector from particles entering from the bottom or 
sides, (bottom) — Detector grid connection for the xenon detector 
HED3 (also HED1). The propane layer rate Is read out as rate V2. 

Figure 4.3 -- Total count rate vs. the total anticoincidence rate for the 
xenon detector. The slight upward trend Indicates that .59£ of the 
particle events are being counted as X-rays. The linear correlation 
coefficient is r. 

Figure 4.4 — The total anticoincidence rate vs. Mcllwaln L. 

Figure 4.5 — Total count rate vs. Mcllwaln L for the xenon detector. An 
upward trend is visible In the data. 

Figure 4.6 — Histogram shows number of times a given value of Mcllwaln L 
occurs. A major frame corresponds to a 40.96 sec readout period. 

The distribution clearly shows a tall to higher values of L (where 
particle contamination Is also higher). 
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Figure 4.7 — The count rate in the second layer of the xenon detector for 
a single pass over the Crab nebula. Detector offsets were 
determined from the sum of many passes. 

Figure 4.8 — The excess noise (above Poisson) plotted as a function of bln 
size for two detectors. The plot Is the sum of 46 different blank 
sky observations. The dashed curve was calculated (IncorrOctly) 
with the weights based on the observed number of counts. For the 
solid curve the calculation was based on the expected number. 
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5. OBSERVATIONS OF ACTIVE GALAXIES 

5.1 Chapter Overview 

In this chapter, adapted from Tennant and Mushotzky 1983, the 
variance analysis of Chapter 2 Is applied to the pointed observations of 
active galaxies. It Is argued that most of the excess variance In the 
active galaxy observations Is due to detector motion over fluctuations In 
the X-ray background. This noise term Is seen In the observations of the 
blank sky. The most Interesting source was NGC 6814 and Chapter 6 Is 
devoted to a discussion of the variability of this source. NGC 3227, NGC 
4151 and MCG-5-23-16 all showed variability consistent with a one day time 
scale. The last section of this chapter discusses the Implications of the 
general lack of variability on time scales of less than one day. 

5.2 Variance 

For all HEAO-1 pointed observations at active galactic nuclei 
Oj/I was determined for bln sizes ranging from 5.12 sec to 1024 x 5.12 
sec m 87 min. The quantity Oj/I Is the standard deviation of the source 
divided by the source flux. This quantity measures an Intrinsic property 
of the source and so is Independent of distance. The results are 
summarized In Table 6.1. The first four columns contain information about 
the observation; column (1) gives the source name In alphabetic order, (2) 
the position, (3) the date of the observation where day 1 = JD 2,443,144.5 
= 1977 January 1, and (4) the length of the observation in hours. Column 5 
gives the mean flux for the observation in m1111-counts/cm 2 -sec. 

Throughout this chapter I report fluxes In these units. A flux of 1 in the 
xenon detector corresponds to 1.0 x 10“*1 erg/cm 2 -sec In the 2-10 keV band 
or 1.8 x 10" 11 ergs/cm^-sec In the 2-20 keV band assuming a power-law 
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TABLE £ I 


(1 ) 

(2) 

(3) 

(4) 

(5) 

(0) 

(7) 

(81 

(9) 

OBJECT 

RA-DEC 

DATE 

LENCTH 

FLUX 

L 

cyi 

q/I 


CEN A 

1322-420 

384.5 

4 

40.5 

0.57 

1.4 

<0.7 

<1.0 

CEN A 

1322-428 

567.7 

It 

36.3 

0.43 

2.8 

1 .8 

3.0 

CYC A 

1057+400 

480.0 

a 

8.5 

95 

5.8 

<3.3 

<0.8 

ESO M1-C55 

1017-580 

466 4 

0 

3.4 

10 

<8. t 

<5.7 

<0.8 

H 1040-506 

1 040—50 1 

620.7 

0 

4.2 

0.28 

11.3 

10.0 

1 .0 

IC 4320A 

1340-301 

573.8 

0 

5.1 

4.1 

9.7 

4.6 

1.1 

MCC-2-S8-22 

2302-000 

523.2 

6 

4.3 

34 

<6.0 

<4.0 

<0.8 

MCC-5-23-10 

0045-307 

407.7 

a 

8.1 

1 .9 

7.0 

5.7 

2.0 

MCC 8-11-11 

0551+464 

643.7 

9 

3.5 

5.1 

9.1 

0.8 

1 .6 

MK M2 

1022+510 

000,4 

s 

1 .6 

11 

<21 .7 

<13.6 

<0.7 

MK 279 

1352+000 

079.4 

13 

1 .2 

3.0 

36.5 

15.7 

0.9 

HK 270 
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7 

1.2 

3.0 
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<1 1 .2 

<0,7 

MK 335 

0003+100 
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13 

1.1 

2.4 

31 .2 

18.3 

8.0 

MK 421 
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513, 1 

e 

1 .8 

0.0 

20.9 

1 1 .6 

0.8 

MK 421 
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7 

1 .2 

4.0 
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<0.0 

I.K 404 

1353+388 
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10 
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13 
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<0,0 

MK 501 

1 052t3SB 
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6 
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10 
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MK 500 
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0 
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20 

5.4 

3.7 

0.8 

MK 500 
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4 

5.3 

23 
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<1 .3 

MK 500 
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657.4 

18 

5.4 

24 

<4.3 

5.3 

1 .3 

MK 500 
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060.3 

10 
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20 

7.3 

8.4 

1 ,4 

MK 500 

2041-100 

666.4 

12 

3.0 

10 

7.0 

6.0 

1 .0 

MK 500 
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6 

1 .1 

2.8 

<19.8 

<1 1 .5 

<0.6 

M 82 

0951+600 
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5 

2.5 

0.007 

<7.4 

<3.0 

<0.S 

NGC 526A 

0121-353 
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0 

2.1 

2.4 
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<6.2 

<0.8 

NCC 931 

0225+31 1 

585.0 

6 

1 .5 

1 .4 
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<0.7 

NCC 2110 

0548-075 
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5 

5.0 

0.87 

<5.3 

<3.3 

<0.8 

NCC 2002 

0043-141 

500.4 

3 

7.8 

1.4 

7,4 

<2.5 

<0.9 

NCC 3227 

1020+201 

602.0 

a 

1 .7 

0.064 

25.5 

20.3 

1 .0 

NGC 3783 

1 1 36-375 

370.0 

3 

3,5 

1.0 

<8.4 

<5.5 

<0.9 

NCC 3783 

1 1 36-375 

554. 8 

13 

3.1 

0.89 

14.7 

8.7 

1 .2 

NCC 4151 

1 208+307 

340.5 

3 

10. 1 

0.38 

4.7 

3.4 

1 .5 

NCC 4151 

1208+397 

524.1 

3 

19.3 

0.73 

5.8 

6.5 

5.7 

NCC 4151 

1208+307 

532.1 

4 

22.8 

0.86 

<4.6 

<4.2 

<2.3 

NCC 4151 

1 208+397 

533.0 

3 

18,8 

0.71 

<2.8 

<2.1 

<1 .8 

NGC 5500 

1410-030 

574.0 

6 

5.0 

0.04 

6.0 

S. 1 

1 .5 

NGC 5548 

1415+254 

566.9 

Q 

5.7 

5.4 

<4.9 

<3.1 

<0.8 

NCC 5548 

1415+254 

730.0 

6 

4.3 

4.1 

<5.4 

<3.8 

<0.8 

NGC 0814 

1940-104 

483.0 

7 

4.1 

0.40 

43.6 

43.3 

8.4 

NGC 7213 

2206-474 

511.4 

0 

4.2 

0.49 

<8.4 

<5.8 

<0.8 

NCC 7460 

2300+086 

547.0 

12 

1 .8 

1 .7 

10.4 

13.9 

l .3 

NGC 7582 

2315-426 

518.3 

6 

3. 1 

0.20 

9.2 

9.9 

1 .5 

PKS 0548-322 

0540-323 

638.0 

8 

2.6 

43 

<10.1 

<8.5 

0.7 

PKS 2151-304 

2151-306 

677.7 

5 

8.1 

810 

<5.2 

5.1 

1 .6 

2A 1210+30S 

1210+305 

516.2 

6 

3.0 

180 

<8.5 

<5.4 

<0.B 

3C 111 

0415+379 

6! 1 .2 

0 

3.9 

29 
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<6.4 

<1 .2 

3C 120 

0430+052 

615.6 

0 

3. 1 

12 

<7.2 

<4.7 

<0.7 

3C 273 

1226+023 

533,9 

6 

9.4 

800 

4.8 

3.3 

1 .5 

3C 273 

1226+023 

546. 4 

6 

8.5 

720 

<3.3 

<2.3 

<0.9 

3C 273 

1226+023 

549.3 

8 

9.0 

760 

<3.4 

<2.4 

<1 .0 

3C 273 

1226+023 

551 .2 

10 

8.7 

740 

<3.3 

<2.4 

<1 .0 

3C 382 

1 833+327 

66) .S 

g 

3.2 

38 

12.2 

<0. 1 

<0.9 

3C 380.3 

1845+797 

723.2 

10 

2.2 

25 

<14.1 

<9.4 
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NOTES i 

(1) Alphabetic by source name 

(2) RA (hours and mins) Deo (degrees and tenths! 

(3) Beginning of rbeervat Ion > day of year 1 077 

(4) Length of observation In hours 

(5) Cf e/om^-ksec 

(6) 2-20 Kev luminosity (H-75 km/sec/Mpc) aeoumlng no absorptlonj unlts-10 srgs/seo 

(7) 328s bln size 
(S) 86m bln size 
(0) 86m bln size 
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spectrum with energy Index o = 0.7. Columns (7) and (8) give the values 
of ffj/I, In percent, for bln sizes of ~ 5 min and ~ 90 min. Finally In 
order to compare with the blank sky observations aj/B In percent Is listed 
In column (9). 

In Figure 5,2 Oj/I, obtained using a 328 sec bln size, Is plotted 
versus source Intensity, In Figure 5.1. Since <?j/I Is computed using the 
entire observation, Figure 5.1 Is relevant for time scales ranging from 300 
sec to ~ 3 hr. The solid line In Figure 5.1 corresponds to a constant 
standard deviation of 0.4 m counts/ cm-sec, which was estimated In the last 
section to be the upper limit for systematic errors. Thus any positive 
detections below this line could be due to residual sky plus detector 
noise. Only three objects, NGC 6814, NGC 4151 (one time), and Cen A (one 
time), lie significantly above the curve. 

We have examined Oj/I down to a bln size of 5.12 sec (a detector 
readout time). Because of counting statistics our upper limits at 5.12 sec 
are ~ 2 times larger than the corresponding values at 328 sec. Since there 
are fewer positive detections of variability at 5.12 sec than at 328 sec, 
we conclude that these objects do not show a large source of variance with 
time scales of less than 328 sec. We conclude that large-amplitude, Oj/I > 
10%, short-term variations on time scales 5 < t < 10^ sec are not, a 
characteristic of the X-ray emission from active galaxies. 

As one goes to longer bln sizes, most of the residual 
non-source-related noise Is averaged Into a single bln. Since the 
spacecraft look position generally samples the entire area available to It 
during one spacecraft orbit, sky noise Is greatly reduced for a 90 minute 
bin size. For a 90 minute bln size one also averages over the particle 
background for one spacecraft orbit. A plot of Oj/I versus source flux for 
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a ~ 90 minute bln size Is shown In Figure 5.2. 

We have attempted to determine the systematics for the 90 minute bin 

size by performing the equivalent test for variability during observations 

of blank sky. Of 49 observations of blank sky, 17 were not constant at the 

90% confidence level. Since for all active galaxies, except Cen A, the 

background gives more counts than the source, we would expect a similar 

fraction of the active galaxies to show variability. Thus we expect 18 i 4 

positive detections, when In fact the observed value was 21. Thus, to 

first order, the majority of our positive detections only Indicate that 

there Is a small amount of excess variance beyond counting statistics, not 
2 

Included in our x test for variability. 

If the same excess variance accounts for the Increased number of 
positive detections both for the blank sky and the active galaxies, then 
the magnitude of the variance should help us decide If any detected 
variability is real or not. In Figure 5.3 (bottom) we have plotted a 
histogram of the standard deviation a B In the background count rate B 
normalized by dividing by B for a ~ 90 minute bln size. The white areas 
represent upper limits, and the shaded boxes correspond to observed 
values. We have computed the excess variance for the entire set of 
observations by adding together the excess variance from each 
observation. For the 90 minute bin size this gave o B /B of 0.7%, which Is 
what one would expect If most of the excess noise was due to background 
fluctuations only. Thus the observed distribution is consistent with the 
expected excess variance being due to "confusion noise." 

For observations of active galaxies we expect a slightly different 
distribution. First, when the detector Is pointed at an active galaxy, the 
count rate Is higher (by definition). This means that our sensitivity to 
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non-source-related variance is smaller during observations of active ti 

galaxies. This Is shown by the fact that the distribution of upper limits 

has moved to the right {Fig. 5.3, top). Since the flux from active 

galaxies was first area corrected for effective exposure, to convert 

to Oj/B we multiplied Oj/I by S/B, where S Is the source count rate (source 

flux times mean area). This means we have effectively Ignored corrections 

due to changing area as the detector passes over the source. 

Based on the results of Figure 5.3, we divide the sources Into three 

groups: (1) those with oj/B > 2%, (2) those with 1.5% < a^/B < 2%, and (3) 

those with Pj/B < 1.5%. Sky fluctuations cannot account for the 
variability seen for the three sources in group 2, whereas any variability 
seen in group 3 can be totally due to sky noise. Sources in group 2 are 
suspect and need further checks. In the next section we will consider the 
sources In groups 1 and 2 In greater detail. 

5.3 Light Curves 

The most variable source was NGC 6814 (Fig. 5.4). This source Is 

highly variable and shows a factor of 2 changes In flux on all time scales 

down to a few minutes (see Chapter 6). Table 5.1 shows that pj/I Is near 
43% for all the bln sizes considered up to 90 minutes. Using the published 
light curve from the Einstein Observatory Imaging Proportional Counter ! 

(IPC) observation of MGC 6814 (Tananbaum 1980), we find that Oj/I ~ 21% for 
a 3 hr bln size. This could Indicate either that the dominant source of 

i, 

variance has a time scale <>< hours or that amplitude of variability j ; v 

decreased, perhaps because of a decrease In flux between the HEAO-1 A2 and 

the IPC observations. j 

In Figure 5.5 we Illustrate the difficulties in determining the time ( 

scale for a poorly sampled observation. The light curve is for our most • 

1 

I * 
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variable observation of NGC 4151. In the top panel we fit all six points 

p 

to a straight line and calculate a x‘ of 1 1.12 for 4 degrees of freedom 

p 

(dof), which is an unacceptable fit to the data. (This compares with x * 
39.80 for 5 dof calculated assuming the source was constant.) In the 

p 

second panel we have tried a model of a step function. For this case x * 

5.39 for 4 dof, which is acceptable. Although It appears (to the eye at 

least) that the transition occurred over the three bins nearest our 

location of the step, we point out that a small systematic deviation 

lasting 45 minutes could easily be introduced by a spacecraft orbital 

effect. In the third panel we have again fitted a straight line to the 

data but in this case have thrown out the first point. For the third 
2 

case x * 4.76 for 3 dof. If the source continued to brighten at the rate 
Indicated In the third panel, then it would double Its intensity In 12 
hr. Thus statistically we cannot distinguish between a 12% flux increase 
on a time scale of 15 minutes or the start of a 12 hr flare. We prefer the 
latter interpretation since it is consistent with previous observations of 
variability from NGC 4151 (Mushotzky, Holt, and Serlemltsos 1978; Lawrence 
1980). This ambiguity as to the correct model Is not resolved as one goes 
to shorter bln sizes, for this case at least, since the smaller bin size 
generates no new Information. 

The third most variable source is Cen A (Fig. 5.4). We note that 
HEAO-1 scanning data for Cen A show It to be slightly extended, presumably 
due to weak nearby sources. Marshall and Clark (1981) have reported a 
nearby source which will appear In our field of view. This source 
confusion can explain most of the observed variability. Notice that there 
Is no Indication of a linear change over the 12 hr observation. 

We will now consider the sources with only weak evidence for 
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Intrinsic changes. First, consider the light curve for H1649-595 In Figure 
5.4 (H1649-595 was tentatively Identified with NGC 6221 by Marshall et al. 
1979), This light curve clearly shows that most of the excess variance Is 
due to short-duration "flares" near the end of the observation. For this 
source the three high points are due to a bright confusing source. Figure 
5.6, a contour plot of the X-ray sky around H1649-595, clearly shows the 
confusing source. The source "flared" only when the detector was looking 
at the region of the sky near this source. We have Included H1649-595 In 
our sample as a warning that "statistically real" variability could have 
many causes. 

Although H1649-595 had a true source of excess variance, it was only 
weakly detected as variable by the Oj/l test. This was due to most of the 
variability being due to flares and, as was mentioned In chapter 2, looking 
at the excess variance Is not the most sensitive method to search for 
single bln flares. Since all sources were examined by eye with a 300 sec 
bin size, I can report that there was no evidence for a single bln flares 
on that time scale. Let us now consider shorter bln sizes. A typical 
observation consisted of ~ 1000, 5.12 sec Intervals. Table 2.1 shows that 
a single bln flare would have had to exceed 7.6a, where a Is the 
statistical uncertainty, before It would have been detected more than 50% 
of the time. If the flare size was less than 7.6o then there is a good 
chance we would have missed the event. For a typical source rate of 3 
ct/sec and a background of 14 ct/sec, a 7.6a fluctuation would correspond 
to a factor of 4.6 Increase In the source flux. Simple Gaussian statistics 
show that 90% of the time all 1000 bins would have statistical fluctuations 
of 4.0a or less. For our typical source a 4.0a Increase corresponds to a 
factor of 2.4 rate Increase. 
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Now consider what happens If the bln size Is doubled. A 7.6o excess 
In a 5.12 sec bln would only correspond to a 5.4o excess In a 10.24 sec 
bln. The corresponding value of X, given by (2.4.1), would drop from 1.3 
to 0.92. Since both numbers were calculated for the same observation, we 
find that the Interesting parameter Is the change In X (or the excess 
variance) as shorter bln sizes are considered. Since I have examined the X 
value for bln sizes down to 5.12 sec I can report that there Is nd strong 
evidence for an Increase at short bln sizes. Thus rare 7. 60 events may 
have been missed by the excess variance method; however, the lack of any 
Increase In the excess variance as shorter bln sizes are considered 
indicate that nothing Interesting was missed. 

Both NGC 3227 and MCG-5-23-16 showed evidence for a linear Increase 
In flux. If such a rate Increase continued, then the source flux would 
double In 10 hr for NGC 3227 and 28 hr for MCG-5-23-16. 

Finally, we come to the MCG 8-11-11 observation shown In the center 
panel of Figure 5.4. In this figure we see small peaks near 2.1, 3.9, and 
5.5 hr. Since these points are ~ 1.5 hr apart, l.e., one spacecraft orbit, 
we suspect that these are not real events. A aj/B of ~ 1.6* confirms that 
the variability seen can be due to confusion noise. If one Ignores the 
short-term variability, then a weak linear Increase Is seen which will 
double the source flux In a few days. We are unable to judge the reality 
of such a trend. 

It is interesting to note that all the variable objects in our 
sample, NGC 3227, NGC 4151, NGC 6814, and MCG-5-23-16, are low-luminosity 
objects with L x < 3 x 10^ 3 erg/sec. This agrees with the HEAO-1 longer 
time scale data which show that lower luminosity sources have a greater 
probability of being variable. In addition, with the exception of 
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NGC 6814, the observed variability is consistent with a time scale of ~ 1 
day. This result agrees with the Ariel 5 observation (Marshall, Warwick, 
and Pounds 1981) that many sources show a 1 day time scale. 

Variability on longer time scales will be reported by Mushotzky 
(1983). However, we can compare the count rates for different pointed 
observations of the same source to obtain some information about longer 
time scales. We find that Cen A and NGC 4151 show large changes over 6 
months. Mrk 509, observed four times in 2 weeks, shows a large-amplitude 
change. This was reported by Dower et al. (1980) using the HEAO-1 A3 
data. Our data for days 657-666 show that the time required for AL = <L> 
is 22 days. Thus we confirm the general nature of the variability reported 
by Dower et al. but indicate a slightly longer time scale. 

5.4 Implications of Rapid Non-Variability 

The previous results indicate that X-ray emission from active 
galaxies rarely varies on a time scale of less than one day. In this 
section I will consider the implications of this non-variability. The next 
chapter Is devoted to a discussion of the rapid variability of NGC 6814 and 
the implications of that observation. Finally, In Chapter 7, I try to 
record le NGC 6814 with the general lack of variability. 

A possible reason for a lack of variability Is for the source of the 
X-ray emission to be very stable. In X-ray astronomy stability appears to 
be an exception. Galactic X-ray sources (associated with accretion) are 
known to be highly variable. However, active galaxies are variable on 
longer time scales and of course NGC 6814 varies on sho"* time scales. 
Therefore, if the lack of variability is due to a stable accretion flow 
then one is forced to explain the reason for this stability. The problem 
is not solved but forced one level deeper. 
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The alternate explanation Is to argue that the shortest time scale 
observed is directly related to light travel time across the source 
reylon. This lack of variability on time scales of less than one day 
Indicates that the X-ray plasma Is about 1 lt-day across. In the black 
hole accretion picture there are two ways to create a large X-ray region. 
Either the black hole Is extremely large or the X-rays are produced far 
from the Schwarzschlld radius. 

If the X-rays come from a region only ~ 10 gravitational radii from 
the central object and If the general lack of variability on time scales of 
less than 1 day tells us anything about the size of the region, then the 
central object must have a mass of 10 9 M 0 or greater. Pounds (1979) 
speculated, based on the apparent ubiquity of variability with a 1 day time 
scale, that all active galaxies contained a 10 9 M 0 central object. The 
Eddington limit for such an object Is 10 47 ergs/sec, which is much greater 
than the luminosity of any object In our sample. Of course, if most active 
galaxies contain a 10 9 M 0 object, then one Is forced to explain why NGC 
68K Is so small. It Is possible that most active galaxies contain dead 
quasars (Lynden-Bell 1969), whereas objects like NGC 6814 might never have 
gone through a “quasar-1 Ike" phase. 

If the X-ray plasma Is not confined to the central object, then it 
could fill a large volume of space. In this picture electrons are heated 
via some unknown mechanism, perhaps in a small volume, and then proceed to 
fill a volume of space about 1 1 across. This may occur In a 
two- temperature disk model (Sha f »' ghtman, and Eardley 1976). Since 
the virial temperature of the protons Is much higher than the corresponding 
electron temperature, an accretion disk can have two temperatures. If the 
electrons are not effectively cooled, then Interactions with protons can 
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heat them to temperatures greatly exceeding their vl rial temperatures. 
Llghtman (1982), considering the possibility of relativistic electron 
plasmas, listed several nongravltatlonal ways to confine the electrons. 
Since It Is notoriously difficult to confine plasmas via electromagnetic 
forces, It Is entirely possible that the central object Is boiling off some 
matter. Thus the X-ray emission can come from two components, the region 
near the black hole and the extended volume. 
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FIGURE CAPTIONS 

Figure 5.1 — The excess noise, Oj/I, Is plotted vs. source flux for ~ 5 
minute bln size. The curved line Is an estimate of the noise 
Introduced by source confusion (see text). The solid dots are 
positive detections of source variance and the bars represent upper 
limits. 

Figure 5.2 — Same as Figure 5.1 except now the bins are about ora 
spacecraft orbit long. 

Figure 5.3 (top) — A histogram of Oj divided by the background rate for 
our observations of active galaxies. The shaded regions represent 
positive detections of variability, whereas the white areas above 
the shaded regions represent the 90% confidence upper limit for the 
nonvariable objects, (bottom) — Same as top except In this case no 
HEA0-1 sources were In the field of view. 

Figure 5.4 — A sample of nine of our X-ray light curves. The top six 
represent the more variable objects In our sample, whereas the 
bottom three represent typical light curves. In all cases a 20 
minute bin size was used to construct the light curves. In the 
upper left of each panel the first number Is aj/I, and the number In 
parentheses Is a^B. Both Oj/I and aj/B are constructed for the bln 
size in the plot. The best fitting linear trend Is indicated by a 
solid line for MCG-5-23-15 and NGC 3227. 
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Figure 5.5 — The most variable observation of NGC 4151 has been fitted to 
three models: (top) a straight line through all the points, 

(middle) a step function, and (bottom) a straight line through the 
last five points. 

Figure 5.6 — The region of X-ray sky around H1649-595. (top) Contours are 
every 10 ct/sec with the lowest level at 10 ct/sec. (bottbm) 
Contours are plotted every .5 ct/sec with the highest level at 10 
ct/sec. For both figures the data has been smoothed to reduce the 
statistical noise. H1649-595 only appears In the bottom plot and is 
at the center of the panel. 
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6. RAPID VARIABILITY OF NGC 6814 

6.1 Chapter Overview 

Material In this chapter first appeared In Tennant et al. 1981. 

First It Is established that the X-ray flux from NGC 6814 varies by a 
factor of 2.5 on a timescale of 90 min. This result by Itself Is one of 
the strongest pieces of evidence for short term extragalactic variability 
reported. The fact that the mean for an entire orbit varied indicates that 
It Is unlikely that near Earth effects could cause the signal. By 
examining shorter time scales we see evidence for large changes on time 
scales as short as 100 sec. Extensive tests for spectral changes showed no 
obvious Indications of a change In the shape of the spectrum. Finally, the 
rapid variability observed from NGC 6814 Is used to constrain various 
models for the X-ray emission from this source. 

6.2 Variability on Long Time Scales 

Examination of Figure 6.1 shows that the flux from NGC 6814 Is 

strongly variable. On an orbit-by-orbit (1.5 hr) basis, the flux varies by 

2 

a factor of ~ 2.5 (Table 6.1). In Table 6.1 we also show the x calculated 
on the assumption that the source was constant during each Individual 
orbit. For all but one of these orbits, the value of x 2 allows us to 
reject the hypothesis of source constancy on time scales less than 1.5 hr 
at > 99.99% confidence. This Is clearly shown In Figures 6.2 and 6.3 In 
which a prominent flux changes on a time scales of less than 500 sec are 
evident. 

Having established source variability on time scales less than 1.5 
hr, we now attempt to characterize the nature and time scale of this 
variability. 
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TABLE 6.1 
NGC 6814 


Orbit 

Mean Flux 
mcts/cm 2 -sec 

2 

X 

(82sec bins) 

Degrees of 
Freedom 

1 

5.6 

279.6 

21 

2 

2.1 

67.9 

23 

3 

4,8 

77.2 

11 

4 

2.5 

12.3 

9 

5 

2.9 

35.9 

6 


6,3 Short Time Scales 

The auto-correlation function (ACF) was computed for the first and 
second orbits using a 10.24 sec bln size. In Section 3.4 we discuss how 
the ACF was generated and how errors in the reported time scales were 
determined. The Inferred characteristic time scale of the source Intensity 
and the errors are given in Table 6.2 for these orbits and the entire data 
stream. All the data are consistent with a characteristic time of ~ 100 
sec. There were Insufficient data for orbits 3-5 for a meaningful 
auto-correlation function to be calculated separately. In Figure 6,4a we 
show the ACF for the complete data set (orbits 1-5 Inclusive). The 
characteristic time for variability is the slope of the natural logarithm 
of the ACF (Figure 6.4b). We detrend the ACF by subtracting the orbital 
mean from each orbit's worth of data. The detrended ACF function and its 
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natural logarithm are shown In Figures 6.4c and 6.4d, respectively. The 
time scale of 200 sec seen when the data are "detrended" confirms that this 
Is a characteristic time for the entire data set. 

TABLE 6.2 


AUTO-CORRELATION FUNCTION RESULTS 


Orbits Used 


Detrend 


Time Scale (sec) 


1 


No 



2 No 230^ 3Q 
1-5 No 248^ 10 
1-5 Yes 101 !26 


6.4 Limits On Spectral Variability 

We have searched for spectral variability In four different ways: 
(1) We have computed the mean spectrum on an orblt-by-orblt basis. (2) We 
have sorted the data on the value of their flux and have computed a mean 
"hardness" ratio for each value of the flux. (3) We have computed the 
cross-correlation function between two energy bands. (4) We have looked 
for spectral variability during a given event. Each of these methods Is 
most sensitive to different types of spectral change as discussed below. 
6.4.1 Orblt-by-Orblt Variability 
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The best fit spuctra on an orblt-by-orblt basis assuming a power law 
model with energy Index a, norm A x TO" 3 counts/sec, and equivalent column 
density of hydrogen in the line of sight are shown In Table 6.3 (all 
errors are 1 0 ). These values are to be compared with the flux weighted 

, ,+ 2.7 nn 

average for the entire observation a « 0.73 ± 0.22, ^ -ii.o x 10^ cm” 

2 (Mushotzky et al. 1980) (90% confidence error). The value of x computed 
on the assumption that a and N^ for each orbit Is the same as the flux 

weighted mean Is 6.8 for ten parameters of Interest. This Indicates that, 
despite a factor of 2 change In Intensity on a time scale of 1.5 hr, the 
time averaged spectrum did not vary significantly. 


TABLE 6.3 

NGC 6814 SPECTRAL FITS 


Orbit 

A 

a 

N h (x 10 22 ) 

2 

X 

1 

20.4 

n C 7 + 0.23 
°- 67 - 0.21 

5 4+ 3 * 2 
0 - 2.2 

1.5 

0 

c 

7.2 

°* - 0.35 

0 c + 3.8 
2 - S “ 2.5 

4.8 

3 

7.7 

0 35 + ^* 23 
U * - 0.21 

n c + 2.5 
2 - 2.3 

14.0 

4 

5.7 

r, cc +0 . 5 1 

0 , 55 “ 0.07 

q+3.6 

7.7 

5 

6.0 

0 , 58 - 0.10 

q+3 .8 

5.3 


Note— 5 degrees of freedom 


C - 
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If we assume that did not vary during this observation 
(consistent with the values In Table 6.3) and fix at the value derived 
for the entire observation, the formal error In a decreases. In Table 6.4 
we show the best fit values of a for each orbit under this assumption. 

There Is a slight Indication that the source spectrum Is marginally flatter 
when the flux Is higher. In Figure 6.5 we plot a versus the 2-20'keV flux 
on an orblt-by-orblt basis and find that a linear relation of the form a = 
mF + b (where F Is the flux) gives an acceptable fit with m = 

-0.054_q!q 22* trend 1s si 9n1f leant at only 2.5 a, however. 
Alternatively, we place a 3 a upper bound of Aa < 0.37 on an orblt-by-orblt 


basis. 


Orbit 


TABLE 6.4 

NGC 6814 SPECTRAL FITS 


5 +0.07 

-0.07 


,+0.19 

-0.16 


5 + 0. 10 
-0.09 


-+0.24 

- 0.21 


20,4 


Note - Absorption fixed at 4.3 x 10 22 


1 q7 + 0»27 
1 , U / - 0.22 
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6.4.2 Limits on Spectral -Intensity Correlation 

We have accumulated data at E > 6 keV (E^g^) and E < 6 keV (Ej_ ow ) 
on an 82 sec basis. We then sorted the data Into 10 Intensity bins or, the 
basis of E L(W . In Figure 6.6 we plot the average of whenever El ow 
was within the prechosen bln limits, along with lines which correspond to 
power law spectra of slopes a * 2.0, 1.75, and 1.50 with N H fixed at 4.3 x 
10*2 cm" 2 . The data are all consistent with the a B 1./5 line. We note 
that each point In Figure 6.6 represents an Independent measure of the 
power law Index. We conclude that there Is no Intensity related spectral 
variability with an upper limit on a change In spectral Index of Aa < 

0.30. (This method of analysis could hide possible variability If the same 
Intensity state has different spectral Indices and If these values scatter 
about a given mean value which Is also characteristic of the total time 
averaged spectrum.) 

6.4.3 Spectral Cross-Correlations 

A defect of the prior two methods Is the relatively large bln sizes 
(30 minutes and 82 seconds), necessitated by the relatively low flux, 
compared to the characteristic variability time of ~ 100 sec. The 
cross-correlation function allows one to test the entire data stream on 
shorter time scales to see whether the softer photons systematically lead 
or lag the hard ones. We report here the result of two cross-correlation 
analyses, one designed to minimize the detectable time scale, and the other 
to maximize the difference between the selected energy bands. 

Cross correlation of the total flux from the argon detector, MED 
(which Is sensitive to photons with a mean energy of F 3 5.46 keV for the 
NGC 6814 spectrum) with the total flux from the xenon detector, HED 3, (F - 
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8.67 keV) shows that the MED leads HED 3 by 6.3 ± 3.3 sec for a bln size of 
15.3 sec. We also cross-correlated photons of energy < 6 keV (F * 5.17 
keV) with photons of energy > 6 keV (F * 12.50 keV) In the HED3 detector. 
For this analysis we used a bln size of 41 sec, since this made the errors 
due to Poisson noise similar to those In the MED versus HED 3 analysis. We 
found that the < 6 keV photons lead the harder photons by -4.1 ± 9.6 sec. 
6.4.4 Spectral Variability of a Resolved Event 

The flux change shown In Figure 6.2 has sufficient statistics to 
examine spectral varlab'llty within a given event. As can be seen In the 
middle panel, the effective spectral Index a did not change (Aa < 0.30) 
during this event despite the fivefold change In source flux. These data 
also rule out variation In X-ray column as the origin of the short term 
variability. Due to the low flux, the hardness ratios are binned in 82 sec 
bins, and thus we cannot comment on spectral variability on the shorter 
time scale Indicated in the previous section. We conclude that, for this 
one event, spectral variability does not occur on the characteristic time 
of intensity variability. 

6.5 Implications of Rapid Variability 

The rapid variability observed In N6C 6814 can place strong 
constraints on the physical processes producing the X-ray flux and the 
environment In which it Is produced. With respect to the N6C 6814 data, we 
shall discuss general constraints and limits we can place on "standard" 
X-rey emission mr.chanlsms. 

6.5.1 General Considerations 

The average luminosity of NGC 6814 Is 5 x 10^ 2 h" 2 ergs/sec In the 
2-20 keV band, here h Is Hubble's constant measured In units of 75 
km/sec-Mpc. Since the amplitude of the variability Is quite large, we will 
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assume that the change In luminosity Is greater than 5 x 10 42 ergs/sec. If 
the X-ray spectrum extends out to ~ 500 keV, as does the X-ray emission 
from Cen A (Baity et al. 1981), then the average luminosity Is ~ 2.1 x 10 4 ^ 
h" 2 ergs/sec. 

If NGC 6814 Is powered by accretion onto a single compact object, 
then we can set a firm upper limit on the mass of this object by requiring 
that the fluctuation time be longer than the light travel time across the 
Schwarzschlld radius. This gives 


M 


3. 7 

< ml x 10' 


(6.5.1) 


for At = 100 sec. On the other hand, we find that 


M > 9 x 10 4 M 0 /h 2 


(6.5.2) 


by requiring that the luminosity not exceed the Eddington limit. Without a 
detailed model, we cannot restrict the mass range any further. As we will 
see, though, a 10® M 0 object Is consistent with thermal Compton models. 

Using the Fabian and Rees (1979) relation given In section 1.4 
with aL > 5 x 10 42 ergs/sec and At^ ~ 100 sec, we find that n > 3% h“ 2 . 
However, If the variations extend out to 500 keV, then n ~ 13% h“ 2 . If one 
requires n < 10%, then h >1.1. Conversely, If h < 1 then one must consider 
models that allow high efficiencies (Thorne 1974) or argue that one of the 
assumptions used In deriving the Fabian and Rees (1979) relation Is 
violated. This could Imply directed motion of either the matter or the 
radiation. 

If the assumptions In the Fabian and Rees relation are In effect, we 
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can deduce two properties of the source. One condition for a minimum time 
scale Is that the optical depth {due to electron scattering) Is near 
unity. This condition tells us that a thermal -Compton model (see section 
6.5.4) may work. A second condition Is that the electrons rapidly radiate 
all their energy. Since NGC 6814 can stay at a high level for a 
long (~ 1000 sec) time, the electron population must be resupplied with 
energy. When matter falls Into a gravitational potential protons’ gain more 
energy than the electrons. Thus the electrons can gain energy via 
collisions with the protons. It is also possible for some of the proton's 
energy to go Into an electromagnetic field (via an accretion disk dynamo 
for example). In this case the created electric field would accelerate the 
electrons. No matter how the energy Is supplied it is clear that if the 
supply is sporadic then the X-ray emission will be variable. 

6.5,2 Synchrotron Radiation and Synchrotron Self-Compton Emission 

If the X-ray flux Is produced by relativistic particles In a "steady 
state" type process via synchrotron emission, or synchrotron self-Compton 
(SSC) emission, then one requires that the particles "live" long enough to 
fill the emitting region. 

The lifetime of a relativistic electron against Compton losses 
Is t, 13 x 107 sec/yu, where the energy of the electron Is yin c 2 and u Is 
the energy density of the photon field. For L ~ 5 x 10^ 2 erg/sec and a 
size R ~ cat ~ 100 lt-sec, u ~ 2 x IQ 6 ergs/cm 2 . This gives t l ~ (7 /y) sec 
which is a factor of 10 shorter than the light travel time across the 
region for all values of y. We conclude that synchrotron and SSC models, 
without continual injection and/or reacceleration and/or relativistic bulk 
motion, cannot work. 

Cava! 1 ere and Morrison (1980) have considered a model in which the 
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particles are continually reaccelerated and in which the radiation rate is 
controlled by the acceleration rate. Their model, for the parameters 
measured in this paper, Implies electron densities n e ~ 1 x IQ 9 and B 
fields ~ 20 gauss for the case in which the X-rays are primarily due to 
Compton Interactions and the optical photons are due to synchrotron 
radiation, a rather higher value for n e and a lower value of B than 
Cavaliere and Morrison found for NGC 4151 or 3C 273. 

6.5.3 Thermal Bremsstrahlung and Blackbody Radiation 

The luminosity from an optically thin sphere emitting X-rays due to 
thermal bremsstrahlung at a temperature T, Is 

L ~ 2,4 x 10" 27 T 172 n e 2 V ergs/sec, (6.5.3) 

where V is the volume of the source region and n e Is the electron 
density. For T > 10 8 K (as required by the spectral fits [Mushotzky et al . 
1980]), one requires n e ~ 1 x 10^ cm" 3 for V ~ 1 x 10 33 cm 3 . This gives 
an optical depth to Thompson scattering t ~ n e oJi ~ 200. Therefore our 
initial assumption of optically thin emission cannot be correct, and the 
output spectrum must be optically thick In shape if the X-rays originate In 
a sphere of hot gas. However, the X-ray spectrum does not have this shape; 
therefore, we conclude that the X-ray emission, without Compton 
amplification, cannot be purely thermal in character. 

6.5.4 Thermal -Compton Models 

Although simple accretion with thermal emission can explain the 
source of power, we need a more exotic emission mechanism to be consistent 
with both the spectrum and the rapid variability. Another X-ray source 
that shows rapid time variations In addition to a power law X-ray spectrum 
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is Cyg X-l. We win now consider models originally constructed to explain 
Cyg X-l but scaled up to L » 10 4 ^ ergs/sec. 

The more successful models have involved unsaturated 
Comptonlzatlon. For this mechanism to function the source must consist of 
a cloud of hot electrons, at a temperature of kT e , with an optical depth to 
electron scattering T gs less than a few and a copious source of soft "seed" 
photons. Under these conditions, soft photons entering the cloud can be 
up-scattered, in energy, to X-rays but do not spend enough time In the 
cloud to come Into thermal equilibrium with the electrons. Shapiro, 
Llghtman, and Eardley (1976) and Katz (1976) showed that the resultant 
spectrum Is a power law for energies below kT e with an exponential fall off 
at higher energies. Pozdynakov, Sobol, and Sunyaev (1979) and Takahara 
(1980) found qualitatively similar results In Monte Carlo calculations for 
a wide range of temperatures and optical depths. The location of the 
rollover is somewhat model dependent with the result from Shapiro, 

Llghtman, and Eardley (1976) occurring at the lowest temperatures. 

However, all of the calculated spectra rollover significantly by 3 kT e . 
There is no evidence for either a rollover or a Wien peak In the spectrum 
of NGC 6814 or of any other active galaxy where spectra have been measured 
to higher energies. This Implies that kT e > 15 keV in NGC 6814. Since Cen 
A, the source with the best determined spectrum, has a power law out 
to ~ 500 keV (Baity et al. 1981), It is possible that kT e > 300 keV. 

For unsaturated Compton models, x es and kT e are related by 


T 


es 


,650 keV* 1/2 2 

(_ ) - 7 

e 


(6.5.4) 


(Sunyaev and Titarchuk 1980), for o = 0.7 as measured for NGC 6814. Using 
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the above we compute that x fle * 4.0, 1.9, and 0.8 for electron temperatures 
of 30, 100, and 300 keV. Notice that although kT e ranges over a factor of 
10, the product (x x kT ) only changes by a factor of 2. This will be 

v!) 6 

useful below. 

Iri any mechanism Involving Inverse Compton scattering, the harder 
photons are on the average scattered more times than the softer photons. 
Llghtman, Glacconl, and Tananbaum (1978) pointed out that since the harder 
photons are scattered on the average more times, they spend more time In 
the source and so tend to lag the soft during fluctuations. Llghtman, 
Glacconl, and Tananbaum (1978) suggested that the rise time of a 
flare t m , measured at energy E, should be given by 

t^ » t R A In (E/E s ) (6.5.5) 

where t R = R/c (size of source region), A Is a constant on the order of 
unity, and E s Is the energy of the original soft photons. More detailed 
calculations (Payne 1980; Llghtman and Ryblckl 1979) have shown that 

m e c 2 

A = tkT (3 + 2a) » (6.5.6) 

6 

where x = Na-j-R and a Is the power law index of the persistent source. 

Since t x kT e appears as a product, we can easily estimate A to with In a 
factor of 2. If we use kT e ~ 100 keV, we find A ~ 0.6. 

In (6.5.5), tjj, Is unobservable since It Is measured from the 
(unobserved) Injection time of the original soft photons. We set - (1 + 
f) t 0 , where t 0 is the observed rise time of the fare and f Is model 
dependent but only a slowly varying function of E s . Solving for gives 
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(1 + f)t Q 

t R “ OTri (E/EJ* 

9 


(6.5.7) 


Before we can go on, we must estimate E $ . Llghtman et al. (1978) 
suggested measuring the rise time of a flare In two energy windows and then 
determining the ratio of rise times: 


. *ml . 1" (E V E s> 

** w 


(6.5.8) 


Solving for E s one finds that, for E s < 0.1 keV, E s Is proportional to Ej 
and Eg raised to large powers. Thus a small error (1%) in D can result In 
a factor of 10 uncertainty In E s . In addition, (6.5.8) Is only 
approximate, since It uses ^ Instead of the observable time t 0 . An 
observational problem Is that the statistics usually are not good enough to 
determine two different rise times In two narrow windows. Thus, we are 
unable to determine E s accurately, but estimating this number will not 
result In too large an error In t^. 

Since Seyfert galaxies are known to be strong IR sources (Rieke 
1978), we choose E s = 1 eV. This Is consistent with the hypothesis of the 
soft seed photons originating in a large accretion disk which should be 
rather cool. We define t 0 as the time for the Intensity to go from 0.37 
I max to I max . With log (E/E s ) = 3.8, we find f ~ 2 using the function 
given by Payne (1980). Thus 


t R ~ 60 sec. 


(6.5.9) 


which says that the size of the source region is slightly smaller than the 


94 


variation time. If we assume the peak flux comes from ■*» 5 Schwarzschlld 
radii, then we estimate the mass of the compact object as 


m m 


c^60 sec 
"TOG 


* 


1 X 10 6 M 0 . 


(6.5.10) 


If the geometry is canonical Kerr, then the peak X-ray flux comes from 0.8 
Schwarzschlld radii (Thorne 1974). For this case M ~ 8 x 10® M 0 . 

Although It has been Impossible to measure D and thus check the 
consistency of the unsaturated Compton model, another test can be made. 
Photons at E 2 (> Z\) should always lag photons at Ej_ by an amount 


4t m ‘ t m2 - ’w ‘ At R l V E l K (6 - 5 - n) 

This equation Is Independent of the soft photon energy but Is a function of 
the electron temperature. It Is possible to measure &t m , using the 
cross-correlation function, for the entire data stream, not just during the 
rising portion of an event. Using t R = 60 sec and A a 0.6, one 
predicts At m = 30 and 17 sec for the < 6 keV versus > 6 keV and MED versus 
HED 3 data respectively. Recall from the observations In §6.4.3 that the 
measured times were -4 ± 10 sec and 6 ± 3 sec, respectively. Both of these 
values of the lag time are consistent with the observed lag time being 
one-third or less of the predicted value. One possibility Is that t R Is 
smaller than our estimate of 60 sec (see below). Another point Is that we 
are using rather broad energy windows which will "smear" the effect. Until 
we obtain more photons, allowing a smaller energy window, we cannot make 
any strong statements. However, the fact that the observed MED versus 
HED 3 delay Is In the right direction is encouraging. 
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Both Payne (1980) and Llghtman and Ryblckl (1979) predict large (and 
so far unobserved) spectral changes during an Intensity fluctuation. We 
note a few reasons that we might not see such changes. First, Llghtman and 
Ryblckl (1979) suggested that one may Identify the observed time for 
variability with the time scale for it change In the soft photon flux rather 
than t R . Equation (27) of Llghtman and Ryblckl (1979) Is then 
applicable. This equation predicts n e > 4 x 10 12 cm" 3 which would result 
In tjj < lo sec, consistent with our upper limit. However, If we use the 
determination of at m * 6 ± 3 sec and apply (6.5.11), we find t^ ffl 22 ± 

11 sec. Secondly, It Is possible that the output spectrum calculated by 
Payne (1980) and Lightman and Ryblckl (1979) is Incorrect because they have 
not calculated a cruly self-consistent model. Therefore, we cannot at 
present use our observations to critically constrain such models. Finally, 
it Is possible that unsaturated Comptonlzation Is not the correct model. 

6.6 Source Lifetime 

We can set an upper limit for the object's lifetime at Its current 
luminosity by requiring that It cannot have accreted more material than Its 
present mass at an assumed constant accretion rate. If we assume a 10% 
conversion efficiency of matter Into radiation, consl itent with some of the 
theoretical discussions of accretion onto a black hole (see Rees 1979 for 
discussion of the energetics of accretion onto black holes and see Thorne 
1974 for a discussion of possible values of the efficiency), the accretion 
rate for NGC 6814 must be > 0.002 M 0 /yr to account for the observed 
luminosity. This gives a lifetime 

10 6 M o r 

Ti < - — * 5 x 10 w yr, 

L ~ 0.002 M 0 /yr 
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which Is considerably shorter than the Hubble time. We cannot rule out the 
possibility that M > 5 x 10 6 M 0 and that the efficiency Is near 30% (Thorne 
1974), In which case the source could have maintained Its present 
luminosity for a Hubble time. If the 10% efficiency assumption Is correct, 
one must consider alternate scenarios. Possibilities are that the object 
has a duty cycle of T L /t Hubbl e < 0.03 or that the high luminosity is a 
recent phenomenon. We note that ~ 0,02 of all galaxies with M v < -19 are 
active galaxies (Huchra 1977). The value of the "on" duty cycle Is 
therefore consistent with the hypothesis that all galaxies show a 
Seyfert-Hke nature for ~ 0,02 of their lifetimes. 
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FIGURE CAPTIONS 


Figure 6.1 — The X-ray flux from NGC 6814 as a function of time for the 

first three orbits. Data from both the argon detector and the xerjn 
detector were added together. The effective bandwidth Is 3-20 
keV. The bln size is 82 sec. 

Figure 6.2 (top) — Enlarged section of the NGC 6814 X-ray light curve 
plotted In Figure 6.1. Time zero corresponds to 24 minutes in 
Figure 6.1. For this figure the bln size was 20.40 sec. The data 
point at 500 sec occurred during a data drop out. (middle) — The 
hardness ratio defined as > 6 keV flux divided by the < 6 keV flux 
in the xenon detector during the event shown above. The hardness 
ratio was computed every 82 sec. The straight lines represent the 
hardness ratios expected for the given values of a. (bottom) — 

The counting ‘'ate In the offset xenon detector for the above 
event. The offset Is 6° from the other detectors. For the data 
shown x for a constant source .itodel was 42.06 for 42 degrees of 
freedom. 

Figure 6.3 — Enlargement of another section of Figure 6.1. For this case 
the bln size Is 20.48 sec and the aata is centered near the time 210 
minutes of Figure 6.1. The decline into the major dip took only 80 
sec. 

Figure 6.4 — Auto-correlation function for NGC 6814. (a) The function for 

orbits 1-5. (b) The natural logarithm of the curve In (a), (c) The 


"detrended" auto-correlation function for orbits 1-5. (d) The 

natural logarithm of the data In {c ) . 
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Figure 6.5 — Power law Index vs. flux for orbital averages. The energy 
Index computed for each orbit Is plotted as a function of the 
average flux during that orbit. The absorption was fixed at ° 
4.3 x 10 22 atoms In the line of sight. The point near 9 x 10" 11 
ergs/cm 2 occurred during orbit three when the spectral fit was 
unacceptable and therefore the error shown may be underestimated. 

Figure 6.6 — Average "hard" flux for a specified "soft" flux. The line 
labeled a « 0.75 Is the best fit line to the data. The fact that 
the line does not go through the origin Indicates that the 
background subtracted was slightly In error. If the spectral shape 
does not change during an intensity change then the data should He 
along a straight line. 
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7, RECONCILIATION 

7.1 Chapter Overview 

In Chapter 5 we saw showed that, with one exception, the active 
galaxies observed by HEAO 1 did not show any evidence for variability on 
time scales shorter than one day. Chapter 6 was about the exception, NGC 
6814, which shows factor of two changes on time scales as short as two 
minutes. In this chapter we consider other observations of rapid X-ray 
variability. NGC 4051 appears to be similar to NGC 6814 In several ways. 
These two objects are very low luminosity active galaxies but when compared 
with optical data, they are overlumlnous In X-rays. This could be caused 
by the X-rays being "beamed" at us or by the objects being new. New 
objects should have weak secondary emission due to time delays. Lack of 
variability on short time scales could be due to either stability or to a 
large source size. We consider both possibilities and consider the 
Implications of the rapid variability of NGC 6814 on such models. 

7.2 Comparison With Other Observations Of Variability 

The large amplitude short term variability observed from NGC 6814 
appears to be unique. Previous observations of rapid variability could be 
modelled as a single step function for Can A (Del val lie, Epstein and 
Schnopper 1978), 3C 273 (Tananbaum 1979) and NGC 4051 (Marshall et al. 
1983). The very large flares reported from NGC 4151 (Tananbaum et al. 

1978) remain unconfirmed and are likely In error (see Appendix F), 
Variability from quasars OX 169 (Tananbaum et al. 1979) and QSO 1525+227 
(Matilsky, Shrader and Tananbaum 1982) are more Interesting. However, both 
observations are marred by the extremely low count rate, and would have 
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been improved with a longer run on source. 

From the above we conclude that rapid variability Is very rare. 

Even for an individual object only one observation of several shows 
variability. In addition, the variability observed Is often quite 
different than what was seen from NGC 6814. It 1$ interesting that the 
light curve that most resembles the X-ray light curve of NGC 6814 (In terms 
of amplitude and time scale) Is the Infrared observation of the BL Lac, 

OJ 287 by Wolstencroft, Gilmore, and Williams (1982). They saw variability 
in only one observation out of ~ 10 on OJ 287 and presumably a large number 
of observations of other sources. 

If one considers longer time scales, the situation is quite 
different. In section 1.5 it was noted that Cen A and NGC 4151 both are 
often variable with a time scale of a few days. These two objects, which 
have similar luminosities to NGC 6814, vary on ~ 1000 times longer time 
scale. 

7.3 Observations At Other Wavelengths 

In this section we will consider NGC 6814 and see whether It has any 
unusual properties that might help us understand the observed rapid 
variability. 

Rleke (1978) found that both NGC 6814 and NGC 4051 had low, but not 
unusually low, IR luminosities. When he compared his measurements with 
those from Stein and Weedman (1976), the two objects that differed the most 
were NGC 4051 and NGC 6814, although Rleke points out that Stein and 
Weedman' s value for NGC 6814 was probably in error. The difference may be 
considered as weak evidence for IR variability. More recently, Glass 
(1979), reporting on IR observations of active southern galaxies, notes 
that from an IR point of view NGC 6814 (and NGC 3783) are only marginally 
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Seyfert-1 Ike. 

Concerning optical correlations, Yee (1980) reports on the very 
strong correlation between the luminosity In Hp and the "nonthermal" 
luminosity for quasars and both broad- and narrow-line objects. NGC 6814, 
NGC 4051 and NGC 3227 lie near each other on the correlation but are a 
factor of 30 weaker than the next strongest Seyfert I galaxy In Yee's 
sample. NGC 4151, which has a similar X-ray luminosity as NGC 6814, has 
over 30 times the He luminosity. Lawrence and Elvis (1982) have shown that 
the X-ray flux correlates with various optical and IR parameters for most 
objects. However, compared to the sample as a whole, NGC 6814 and NGC 4051 
are underlumlnous In [0 III], 3.5 pm, and 10 pm flux relative to their 
X-ray luminosities. It Is also Interesting that, in the Lawrence and Elvis 
plots, NGC 6814 did not have an unusually low He flux. 

There are several ways that one can underproduce optical line 
emission. The primary energy source for the optical line emission Is the 
absorption of UV photons. Therefore a lack of line emission could be due 
to either a lack of absorbing matter or a lack of UV photons. The X-ray 
spectrum of NGC 6814 (Mushotzky et al. 1980) shows absorption caused by ~ 4 
x 10 22 atoms/cm 2 along our line of sight. This Is similar to the 
absorption observed In NGC 4151 which would imply that both nuclei have a 
similar amount of matter around them. As for a lack of UV photons, this 
would Imply that an object like NGC 4151 has a source of UV photons (which 
NGC 6814 would lack) above the power law continuum. There Is no evidence 
for an additional source of UV photons In NGC 4151 (Malkan and Sargent 
1982). Also, It Is possible that NGC 6814 has an anomolous Ionization 
parameter. 

Another w Ay to decrease the X-ray flux relative to the line emission 
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Is to assume that the observed X-ray flux Is not typical. If the X-ray 
emitting plasma is moving with velocity p « v/c then the integrated 
luminosity will be increased by a factor of 6^ where 6, the relativistic 
floppier factor, is given by 

6" 1 « y(l - 0COS6) 

where y • 1/7(1 - 0 ^) and 0 Is the angle between the direction of motion and 
the direction along the line of sight. In addition the observed time scale 
for variability will also decrease by a factor of 6. The decrease In the 
true X-ray luminosity, relative to the (Inferred) observed luminosity, by a 
factor of 30 w,)uld require a 6 a 2.3 and of course the true time scale 
would be a factor of 2.3 longer. 

A final method to decrease the line flux is to invoke a time delay 
as the X-ray flux is turned on. X-rays can proceed directly from the 
central source to the Earth. However, In order to produce the line 
emission, the Ionizing photons must first propagate from the central source 
to the clouds before a photon from the lie can be emitted which can 
propagate to the Earth. This results in a longer path length for the line 
emission then for the X-rays. Therefore variability from the lines will 
always appear to lag variability from the ionizing source and appear t? be 
smeared in time and amplitude. 

We note that the probability of seeing 1 object In about 40 In the 
process of turning on depends greatly on the totally unknown X-ray lifetime 
of these objects. Consider the following possibility. After the X-ray 
source turns on e it will take some time t ~ R/c to illuminate the entire 
narrow-line region. This will be the turn-on time. The clouds quickly 
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f"iat up and accelerate to v/c *» 1/100 and proceed to escape from the 
nuclear region. In section 7,3.3 It will be argued that thermal 
reradlatlon from these clouds could affect the observed X-rays. Thus the 
"on" time, which starts when the clouds are heated and ends when the clouds 
escape, will last ~ 100 times longer than the turn on time. Therefore U 
is not unlikely to see 1 object in a sample of ~ 4U in the process of 
turning on. Of course, if this is true, then it implies that an active 
nucleus has gone through many on-off cycles in the history of the galaxy. 

The idea of recurred activity from galaxies Is not new. Oort 
(1977) gives evidence of nuclear activity for our gal any 10® to 10? years 
ago. I will give two examples of models which "predict" this time scale. 
Van Bueren (1978) suggested that tidal disruption would fill the potential 
cusp around the central black hole with gas. As the cusp fills the 
radiation pressure would build up. After 10® - 10^ years the radiation 
pressure would exceed the gravitational pressure. Therefore Van Bueren 
predicts a long period of relative quiet followed oy an explosion. Sanders 
(1981) considered the Interaction of the central black hole with molecular 
clouds In the nuclear region. He calculates that a molecular cloud would 
collide with the hole once every ~ 10^ years. Activity would last ~ 10® 
years and so Sanders predicts a 1% duty cycle. In both these models one 
expects to see X-rays only during the short active phase. 

If one turns on an X-ray source in an originally "normal" nucleus 
then one expects first to detect the effects on material closest to the 
nucleus. Since the broad-Hne region is often less than one pc in size, an 
external observer would see the entire region Illuminated In only a few 
years. However, the narrow-line region, which Is 100-1000 pc across, will 
only be partially illuminated for young objects. Thermal reradlatlon from 
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dust exterior to the nucleus should turn on with the narrow-line region. 

In this picture Seyfert 2 galaxies represent the class which Is turning 
off. Thus a large fraction of the observed differences In Seyfert galaxies 
would represent different snapshots of a single process. However tempting 
such a model Is, we note It does have problems. Dust near the nucleus will 
be at the highest temperature and will be observed to turn on first. This 
Is contrary to observations In that NGC 6814 does have 10 pm emission 
(Rleke 1978), Indicating a cool temperature for the dust, and NGC 1068 (a 
Seyfert 2) has an IR spectrum requiring a dust temperature of ~ 1000 K 
(Jones et al. 1977). 

7.3.1 Stability Related to the Eddington Limit 

Cowie, Ostriker and Stark (1978) and more recently Krolik and London 
(1983) considered the long term stability of accretion. If the luminosity 
Is < 1 % of the Eddington limit stable flows can develop. This would 
explain the observed stability for our objects If they are < .01 L EDD . 

They assumed that the central object was Imbedded In a homogenous gas and 
so the shortest time scale obtained was related to the sound travel time 
across the sonic radius. Therefore they Jld not consider the extremely 
short time scale observed for NGc 6814. This model does predict that NGC 
v '14 should be highly variable on a time scale of - 1 year. This is In 
agreement with Hal pern (1982) who reported that the 2-10 keV flux from NGC 
6814 was down by a factor of 10 one year after the HEAO-1 observation. 

7.3.2 Increased Number Of Shots 

In the Gullbert, Ross, and Fabian (1982) model a cloud Is heated via 
some unknown mechanism and then allowed to cool via Inverse 
Comptonlzatlon, Since the cooling times are very shurt, the observed X-ray 
spectrum Is a time average. Their calculated "averaged" spectrum 1* in 
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good agreement wi*h observations of the X-ray spectra of active galaxies. 
However, If the X-ray source Is powered by discrete events, these events 
should glue rise to low-amplitude varfblllty. To Illustrate we apply the 
shot noise model. In this model X events per time Interval earh rise to 
amplitude h and then exponentially decay away with time scale t. Thus 
using I a xht and oj ■ xh t/2 we construct 

tfj/I - (2 /Xt) 1/2 . 

We set N a Xt, which represents the number of "simultaneous" events taking 
place. If we assume that t is longer than the smallest bin size we used 

p 

(see Sutherland, Welsskopf, and Kahn 1978, Appendix A, to see how aj Is a 
function of bln size) and that X Is greater than one shot per day, then we 
can use Oj/I from Table 5,1. For NGC 6814, N ~ 10, which says that at any 
one time on the average, 10 clouds dominate most of the X-ray flux. A more 
typical value of aj/I near 10% Implies that N ~ 200, and for Cen A, N Is 
greater than 1000. Since Cen A, NGC 4151, and NGC 6814 all have roughly 
the same luminosity, accounting for the lack of rapid variability in Cen A 
and NGC 4151 by Increasing the number of shots only works if the shots 
become much more numerous, and as a result, each event becomes much less 
luminous. 

One way to reconcile the Gullbert, Ross, and Fabian (1982) model 
with the lack of variability is to assume that the heating and cooling are 
taking place in a continuous matter. 

7.3.3 ft. owth Of New Source Of Soft Photons 

In the Inverse COmpton reflection model (Llghtman and Ryblckl 1980), 
low-energy (soft) photons enter a region of energetic electrons. The soft 
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photons Inverse Compton scatter from the electrons to emerge from the cloud 
as X-rays. In this section, we will examine the Inverse Compton reflection 
process and find a constraint on the temperature of the soft-photon 
source. We will make no assumptions about the electron population, l.e., 
the distribution can be either thermal or nonthermal. We will assume that 
the X-ray emitting plasma Is quite small and stable. By stable we mean 
that any variability seen Is due to variations In the soft=photon source 
and not due to changes In the plasma Itself. Tennant et al. (1981), using 
the results of Llghtman and Ryblckl (1979), pointed out that the lack of 
spectral change during the Intensity variations observed In NGC '6814 Is 
consistent with this Interpretation. Thus the lack of rapid variability In 
most sources could be explained by the growth of another "stronger", but 
more constant, source of soft photons. If the X-ray source Is slowly 
heating up Its environment, then thermal reradiation could be the new, 
constant, strong source of soft photons. This is consistent with the 
observed deficient IR flux for the rapidly varying galaxies. 

We will now find the minimum temperature that a thermal source can 
have and still provide enough photons (for Comptonizatlon) to generate the 
observed spectrum. If the soft photons are at a temperature of kT, then 
the observed pow^r law will extend from ~ 3 kT cut to energies determined 
by the temperature of the scattering cloud. For the case of NGC 6814 the 
total number of photons radiated in the Comptonized spectrum Is at least 

>1 A 7 y in-3 A a 2 keV r“l*7 AC 

N c = 4.7 X 10 4 tt d J 3 E dE, 

which corresponds to 2.1 x lO^ 2 photons/sec (kT/1 eV)~0-7 a t a distance d 
of 21 Mpc. The number of blackboc|y photons Impinging onto the X-ray region 
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Is given by 

N bb < 2.7 x 10 49 r^ 0 Q (^-y) photons/sec, 

where r^g Is the radius of the X-ray emitting region In units of 100 
It-sec. N bb Is maximum when the dilution factor Is equal to unity. Since 
for Compton scattering, photon number Is conserved, we set N c * N BB and 
find that 

kT > 6.0 rj^ 2 eV. 

If r^QQ < 1, then we find the thermal source must have kT > 6.0 eV (70,000 
K). Dust grains would quickly vaporize at this temperature. This leads us 
to conclude that thermal radiation from dust cannot provide enough photons 
to generate the observed spectrum for a small X-ray source. This problem 
Is serious for NGC 6814, where the "hot spots" which provide the soft 
photons must be small and few In number In order to account for the rapid, 
large-amplitude variability seen. If these spots reside outside the X-ray 
region, then the dilution factor must be very small, and hence kT » 6.0 
eV. 

7.4.1 Large Source Size 

To summarize the results of the last section. If the source of soft 
photons Is thermal in nature and If some of the optical and/or IR emission 
comes from the X-ray plasma directly, then the X-ray cloud must have r 100 
» 1. Since there appears to be some correlation between 3.5 pm IR 
emission and X-r<v flux for most active galaxies (Lawrence and Elvis 1982) 
and since a large source size Is consistent with the absence of rapid X-ray 
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varf ability reputed here, we will consider the possibility of large X-ray 
regions. In order for 1000 K blackbody photons to be the soft-photon 
source, we find that r^o > 2700 (" 3 lt-days). This size Is consistent 
with previous observations of variability (Marshall, Warwick, and Pounds 
1981). An X-ray plasma this large could be generated In one of two ways In 
the black hole accretion picture. Either the central object Is large, 
hence very massive, or else the X-rays come from a large region not 
directly related to the central object. 

7.4.2 Two Components 

In this section I assume that X-rays are produced In two components; 
a small variable component that would dominate In NGC 6814 and large more 
constant component that would dominate for the remaining objects. Let us 
assume that the total power produced by a galactic nucleus comes out In two 
forms— the Immediate production of X-rays Lj and some Initially unobserved 
power P. The latter could be In the form of relativistic electrons, as 
mentioned above, or in v-rays, as In the Penrose photoproduction model 
(Lelter 1980). The quantity P will slowly fill a reservoir with energetic 
electrons. When steady state Is reached, the luminosity of the reservoir 
will be <P>. Therefore, the total X-ray luminosity L x will be 

L x * Lj + a<P> = L^ + r<L^>, 

where a is the fraction of the reservoir's luminosity which comes out as 
X-rays and r Is a<P>/<L^>. Below we will assume r Is constant and 
that a ~ 1. 

Since only Lj will show rapid variability, let us consider what 
happens to Oj/I when one adds a variable source to a source of constant 
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Intensity I 0 . If we assume that the Intensities of the variable component 
are uniformly distributed between 0 and Al, then 

oj/l = AI[3 1/2 (2I 0 + Al)] 


for the sum of both components. We will now assume that the reservoir does 
not vary on the time scales we have sampled and also that the X-rays 
produced near the central object are highly variable. If we set I 0 * r<L^> 
(the luminosity of the reservoir) and Al * 2<L^> (the range of luminosities 
for Immediate X-ray production), then 


Oj/I 



_1 

(1 + r )* 


We define F to be the fraction of the total X-ray flux coming from the 
compact-variable region, which Is given by 1/(1 + r) = 3^ 2 a^/l * F. 

We can check the consistency of this model by assuming that F = 75% 
and that the average flux is 4 for NGC 6814. Thus, If the assumption of a 
uniform distribution of Intensities Is correct, we would expect to see the 
source vary from ~ 1 to ~ 7. The NGC 6814 flux shown In Figure 6.1 varies 
from ~ 2 to ~ 7, which Is In rough agreement with our model. 

Since our upper limits for a^/I typically lie In the range of 10%, 
we find that F Is typically less than 17%. For NGC 4151, which has 
<jj/I < 6%, we find that F < 10%. Thus the constant component would have to 
grow by a factor of 10 (relative to the variable component) In order 
for Oj/I to decline from the NGC 6814 value to the level observed for 
NGC 4151. One possibility Is that NGC 6814 has not filled Its reservoir 
and 'thus has not come to steady state. This Is unlikely since It Implies 


% 


116 


that when steady state is reached, the luminosity will be 10 times what It 
Is now. As shown in Tennant et al. (1981), such a high total luminosity, 
if it were variable on the same time scale, would clearly violate the 
Fabian and Rees relation (1979). However, it is unclear as to whether the 
Fabian and Rees relation applies in this case since we are talking about a 
steady state condition. If the X-rays from the variable component pass 
through the reservoir, electron scattering could reduce the amplitude of 
variability. If the electrons have the correct power-law distribution, 
scattering will not greatly alter the spectrum. One could also argue that 
NGC 6814 does not have a reservoir for some unknown reason or that the 
X-rey production efficiency a for the reservoir is low. 
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APPENDIX A 
Weighting The Data 

In this section we will consider the effects of basing the weights 
on the observed number of counts rather than the expected number. We will 
assume that the observed count rate has a Poisson distribution with a mean 
of p. In order to consider analytic results, we will assume that u is 
small. Since p Is the mean number of counts per bln, It Is always possible 
(subject to experimental considerations) to generate data with small p by 
binning the data Into a large number of bins. 

First, we will need to determine the appropriate weight for a bln 
that contains no counts. To do this consider dividing the data Into such a 
large number of bins that each bln contains either zero or one count. For 
this case 


- °- w o f o * 1 - w l f l 

w ' Vo + V 1 


(Al) 


where f n Is the fraction of the total number of bins that contain n counts 
and w n Is the corresponding weight. 

We know that, on the average, f n Is the probability of getting n 
counts. Since p Is small, we can expand P p (x; p) to first order to obtain 


f 0 = P p (0; p) = 1 - p 


fl = P p d; u) » p 

f n = Pp{n; p) = 0, If n > 1 


(A2) 
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When we apply (A2) to (Al) we find that 


<x> 


V , V 

w “ W 0 ll-M) + W 1 P w 0 + pYwy w q T 


(A3) 


If we require <x> w = p we find two solutions; either p«l or w 0 = w^. 

Clearly p=l does not apply to our case and so we find w 0 * wj_. Thus, If w n 
= 1/n we find that a good definition for w 0 Is to let w 0 a 1. 

Now assume that we have some bins with 2 counts In them but 
(effectively) no bins with 3 or more counts. Now we must expand P p to 
second order which gives 


Pp(0; p) 
P p (l; p) 
Pp(2; p) 


~ .. . p 

= 1 - p + - 


* y - v 
' T 


2 

T 


The weighted average Is given by 


. . 0 . w 0 f 0 + ! . w 1 f 1 + 2 . w 2 f 2 

w Vo + vr<^ 


where f n Is the fraction of the bins that contain n counts. Again, for a 


large number of bins, f n = P p (n; p), and so 


<X>u, * 


w^(p - p^) + 2 w g (p 2 /2) 


W Wp(l - p + p /2) + w^(p - p ) + w^(p /2) 


Now we make use of the fact that w 0 = w^ = w, and obtain 
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WM + U 2 (W-W 2 ) 

<x> w * H 2 

{1 ** TJ — ) W + 7J— W 2 

It Is clear that w 2 * w gives the correct result of <x> w ■ \i. However, If 
the weights are based on the observed number of counts (Wj * w 0 * 1 and w 2 
s 1/2) then 

* U (1 “ • 

This shows that a bias exists when the weights are based on the observed 
number rather than the expected number of counts. 
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APPENDIX B 

Mcllwaln's L Parameter 

In this appendix I briefly describe what Mcllwaln's L parameter Is 
and wKy It Is useful. 

In a pure dipole magnetic field the trapping of charged particles Is 
relatively simple. Since the particles follow the field lines, a useful 
coordinate system Is one that labels field lines. In such a system, a pure 
dipole field Is described by 




1/2 


R ■ L 


2 

cos A 


(B.l) 


where (R,A) are the polar coordinates (radius and geomagnetic latitude) of 
the point In question, B is the magnetic field strength, and M Is the 
dipole moment. The quantity L Is constant along a field line and measures 
the distance from the origin of the field to the point at which the field 
line crosses the equatorial plane. Figure B.l, adapted from Mcllwaln's 
1963 paper, illustrates this transformation. 

The Earth's field Is not a pure dipole. Mcllwaln (1963) proposed a 
coordinate system which preserves the simplicity of the (L,B) system 
described above. In effect, one maps the Earth's field Into a pure dipole 
while preserving an adiabatic Invariant. In the new coordinate system the 
observed particle flux should be a simple function of (L,B). Figure B.2 
from Mcllwaln's paper, illustrates where the trapped particles lie in the 
(L,B) system. In Figure B.3 I have Illustrated the location of the HEAO-1 
spacecraft In the (L,B) plane for a typical pointed observation. This 
figure shows that HEAO-1 Is below the bulk of the radiation belt. Electron 
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events occur when a group of trapped particles loses some energy In the 
Earth's atmosphere and starts to follow lower field lines. These events 
tend to be of short duration, since electrons are not trapped at these 
lower altitudes. We have found that such events tend to be seen at high 
values of L although Appendix C mentions two events that were seen at lower 
values of L. 

E. Boldt (private communication) has pointed out another effect. 
Primary cosmic rays leave tracks throughout the detector and so are 
eliminated by the anticoincidence logic. However, cosmic rays can produce 
secondaries which will not be eliminated. For example the cosmic ray could 
knock out a K-shell electron in the walls of the detector. When the atom 
recombines It Is possible for It to emit an X-ray Into the detector volume 
which would contribute to the background rate. It Is also possible for the 
cosmic rays to produce y~rays (via neutral pi on production in the 
spacecraft). These y-rays can produce low energy Compton electrons in the 
detector volume. If the electron deposits its energy In one cell then it 
will not trigger the anticoincidence logic. Thus we find that some 
fraction of the detector background should be related to the particle 
background measured by the anticoincidence rate. 

Mcllwaln L Is a good measure of cosmic ray flux (see Smart and Shea 
1967). The reason Is simple to understand. The Earth's magnetic field 
tends to keep charged particles out. Near the magnetic pole (high L) 
cosmic rays from some directions can follow a field line down to the 
Earth's surface. Near the equator (low L) cosmic rays must cross field 
lines no matter what direction they come from. Therefore low energy cosmic 
rays are excluded from these regions of low L. Since the cosmic 
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ray spectrum falls rapidly with energy^ excluding the low energy particles 
results In a total particle flux decrease Thus high values of L 
correspond to higher cosmic ray flux and hence Increased contamination. 


123 


FIGURE CAPTIONS 


Figure B,1 — The mapping of the geomagnetic coordinates R and \ onto the 
B, L plane according to (B.l). 

Figure B.2 — Contours of trapped particle flux In the B, L plane (Mcllwaln 
1963). 

Figure B,3 — The location of HEAO-1 detector Is plotted on the B, L plane 
every 40.96 sec. 
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APPENDIX C 
The Rejects 

In this appendix I discuss the Instances of apparent variability 
which were not actually associated with the target sources. 

When searching for source variability, I examined the rates for all 

2 

mlnus-y points, and used x to test whether the rates were constant. If a 
source was flagged as being In the field of view then any detected 
variability was assumed to be due to the source. No attempt was made to 
determine If the source actually causeu the variability or not. The rates 
were also examined by eye. Two flares were originally discovered using 
this method (see section 4.6). 

The first flare was observed on duly 11, 1978. The light curve Is 
displayed In Figure C.l. The spectrum was unusual In that It rose with 
energy. We were also concerned that the two xenon detectors (HED2 and 
HED3) did not see the same flux. The offset xenon detector (HED1) did not 
see the event at all. When we later discovered that the flare was 
coincident In time with a solar flare, It became obvious that the counts 
were due to low energy y-rays entering the detector from the side. 
Attenuation In the walls of the detector greatly reduced the flux of < 20 
keV flare photons. HED3 was closest to the sun and was the detector with 
the largest count rate. HED1 was behind both HED2 and HED3 and so showed 
the smallest (l.e. no) response. Later In the day a second, larger solar 
flare occurred. The rate from this event Is Illustrated in Figure C.2. 
This gigantic flare took place while the spacecraft was scanning and so Is 
presented here for Its curiosity value. 

The second mlnu/j-y flare was observed on December 29, 1978. Its 
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light curve Is shown In Figure C.3. It is Interesting In that It resembles 
the light curve for the first solar flare. The spectrum, though, Is well 
fit by c simple (decreasing) power law that Is typical of many X-ray 
sources* All three detectors which observed the source saw the same flux 
whereas the off source detector saw nothing. We believe that this event Is 
a "nearby" galactic source (see Tennant and Swank 1983). Models of 
accreting neutron stars predict such events. 

During the observations of active galaxies there were two events 
that we do not consider to be due to a galaxy* These events are 
Illustrated In Figures C.4 and C.5. In both cases the offset detector saw 
the events. This Is a strong Indication that the flares were caused by 
electron contamination. Moreover, as shown In the bottom half of Figures 
C.4 and C.5, the propane veto rate was extremely high. The propane only 
weakly responds to X-rays but does respond to the passage of charged 
particles. This Is also strong evidence that the events were particle 
Induced. Low energy electrons would only have been detected in one (mostly 
the first) layer, and so could have produced these events. 

The reason these events slip through our electron flag Is quite 
simple. We calculate the number of electrons stopping In the detector 
based on the total anticoincidence rate. However, If the electrons are 
predominately stopped in the first layer, then they will not substantially 
increase the anticoincidence rate. Thus, for this case, we underestimate 
the number of electrons entering the detector. Examination of the PHA data 
for the events seen In Figures C.4 and C.5 shows that the contamination Is 

mainly due to a copper fluorescence line. Therefore, the electrons are not 

strongly affecting our count rate directly but rather are knocking out 

K-shell electrons from the copper collimator. X-rays are emitted when the 
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copper atoms recombine. It Is Important to realize that flare-like events 
In our data are very rare. These were the only two events seen in the more 
than 150 hours of data examined. To be on the safe side, it Is recommended 
that future work use an electron flag based on the first layer rate Instead 
of the total anticoincidence rate. 


FIGURE CAPTIONS 


Figure C,1 — Xenon total rate vs. time during a solar flare. Time 0 
corresponds to 15980 sec. 

Figure C.2 — Total rate In the xenon detector vs. time during a major 

solar flare. The tic marks are not related to the size of the error 
bars. 

Figure C.3 — The sum of the xenon and argon detector rates vs. time during 
an X-ray flare. 

Figure C.4 — (top) The rate In the top layer of the xenon detector vs. 

time. An excess Is seen around 77000 sec. (bottom) The rate In the 
propane layer In front of the xenon. 

Figure C.5 -- Same as C.4 except this event occurred at 30000 sec on day 
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APPENDIX D 

Variable Soft X-ray Flux From the Pointed Observation at ESO 141-G85 
In the top half of Figure D.l, I have plotted the argon detector 
count rate for photons with energy < 3 keV. Figure D.l (bottom) Is the 
corresponding plot for energies > 3 keV. It Is Important to realize that 
both plots contain data from exactly the same times and from the same 
detector. This data was from the same observation of ESO 141-G55 In which 
Hushotzky et al. (1980) reported a soft (< 3 keV) excess. Therefore taken 
at face value the figures Indicate that this soft excess Is showing 
variability on a time scale of only a few hours. From an astrophyslcal 
standpoint the luminosity from ESO 141-G55 Is about 10 times greater than 
NGC 6814 and the Inferred time scale Is about 10 times longer. 

The first Impression that one obtains looking at Figure D.l Is that 
the variability Is not due to any near-Earth effects since there Is no 
strong orbital dependent signal In the data. There are several other 
reasons why It Is unlikely that the soft excess Is due to particles. 

1. The particle flux and spectrum would have to be very unusual to give 
rise to counts In only the < 3 keV band. 

2. Since the particles are not vetoed they must deposit most of their 
energy In the first layer. Figure D.2 shows the same soft 
enhancement In the second layer of the argon detector. Again one 
sees the variable signal In the low energy window, while the rate In 
the high energy window Is roughly constant. 

3. The xenon detector has a very weak response to events < 2.5 keV and 
so one cannot clearly see the time signature of the event. However, 
the total spectrum of the event Is consistent with the xenon 
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detector seeing the same soft excess that the argon detector sees. 

4. The xenon detector has a top layer propane veto that responds mostly 
to particles. The rates from this veto laye* did not show the 
characteristic time signal seen In Figure D.l. 

In summary, the signal that we see Is consistent with both detectors 
responding to X-rays. Since the two detectors respond differently to 
particle events we conclude that It Is impossible for particles to give 
rise to the observed signal. 

However the variability Isjiot from ESO 141-G55. During this 
observation the spacecraft spent some time pointing 6° aw ay from the 
source. In Figure D.3 which has the same scale as Figure D.l, data Is only 
included when the source Is not In the detector's field of view. Again one 
should notice that only the < 3 keV rate Is affected. The fact that the 
characteristic time signature Is seen Indicates that the soft X-rays are 
coming from an area extending over 6° on the sky. The fact that the rate 
varies Indicates that the source must be close to the Earth. Thus the 
observation could be explained by some sort of high altitude air-glow. 

In a quest for more data, I regenerated the data file without 
selecting on Mcllwaln L. The resulting light curve Is displayed In Figure 
D.4. The high bins all show a clear orbital signature which one would 
expect If the effect Is produced near the Earth. What Is surprising Is 
that the high rates occur when the source Is coming out of Earth 
occultation but the rates are low going Into Earth occult. Since ESO 
141-G55 was In the morning sky, for the observation In April 1978, we find 
the detector was looking at the sunlit Earth going Into Earth 
occultation. When the source came out of occultation not only was the 
detector looking at the dark Earth but the Sun was below the spacecraft's 
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horizon. 

It Is Interesting that selecting data with low values of Mcllwaln L 
greatly decreased the soft excess. This Is mostly due to chance. Both 
Mcllwaln L and the soft excess are related to the spacecraft's position 
about the Earth. Therefore, If the soft excess occurred at high values of 
L during one orbit, It would continue to appear at high L values during 
later orbits. The large peaks shown In Figure D.4 do not appear In the 
background. When Mushotzky et al. (1980) constructed the spectrum of ESO 
141-G55 they did not select low values of Mcllwaln L. As a result, when 
they subtracted the off-source spectrum, the soft excess was still 
visible. When I reconstructed the on-source minus off-source spectrum 
using only data at low values of L, the soft excess was no longer seen. We 
conclude that there Is no longer any evidence for a soft excess In the 
spectrum of ESO 141-655. 

An explanation tliat fits the data Is that the soft ((-rays are 
produced by a very high altitude air glow. Previous air glows that were 
detected with HEAO-1 were excited by solar UV radiation and so were 
strongest on the sunlit side of the Earth and rarely seen above 100 km. We 
routinely reject data If the lower 200 km of the Earth's atmosphere Is 
anywhere In the detector's field of view. To my knowledge this Is the 
first HEAO-1 observation of air glow that Is both above 200 km and on the 
night side. 

The only thing that I can think of that excites the upper atmosphere 
at night Is trapped particles entering the atmosphere. Thus what we were 
observing was an aurora australis In the X-ray band. Since ESO 141-G55 Is 
located rather far south In the sky, the spacecraft was looking In roughly 
the right direction In order to see an aurora. More Importantly, the 
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geophysical magnetic Indices show that the ESO 141-G55 observation occurred 
on the most disturbed day In April 1978. In fact the Indices were oi.ly 
slightly higher In early May when there were extensive reports of auroral 
activity. 
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FIGURE CAPTIONS 

Figure D.l — (top) The soft (< 3 keV) count rate In the argon detector Is 
plotted vs. time during the observation of ESO 141-G55. An excess 
Is seen with a peaking In the period 2 to 5 hours. For this and the 
remaining figures In this appendix, time 0 corresponds to 30,488 sec 
on day 466 (of 1977). (bottom) The corresponding hard flux for the 
argon detector. The excess Is not visible but some variability, due 
to ESO 141-G55 not being area corrected, is seen. 

Figure D.2 — Same as figure D.l except this time counts from the second 
layer of the detector are plotted. A weak excess Is seen from 2 to 
5 hours In the soft flux (top) but not In the hard (bottom). 

Figure D.3 — Same as figure D.l except now data is plotted when ESO 

141-G55 Is not In the field of view. Again a flare Is seen in the 
soft flux (top) but not the hard (bottom). 

Figure D.4 -- Same as Figure D.l except now data from all values of 

Mcllwaln L are Included. One can now see an up to 6 ct/sec excess 
In the soft flux (top). Tho flares occurred when the detector came 
out of Earth occultatlon. Some of the noise seen In the hard flux 
(bottom) is due to the small Increase In the background when 
Mcllwaln L Is high. 
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APPENDIX E 


The Future of Timing 

The future of X-ray astrophysics requires detectors with larger area 

and reduced Internal background. It Is Interesting that the next two 
8 

major (> 10 $) Instruments optimize one or the other of the two basic 
requirements. The X-ray Timing Explorer (XTE) Is designed to be an 
Inexpensive (~ 10® $) large area experiment. To meet these requirements It 
Is proposed that XTE be built with an ~ 1 m 2 (10 4 cm 2 ) proportional counter 
array. On the other hand, AXAF {Advanced X-ray Astrophysics Facility) Is 
designed to have very high resolution and thus a low background (per 
pixel). Due to the high cost of the telescope, the total collecting area 
Is a more modest 2000 cm 2 . Although both Instruments will be quite useful, 
I believe that the mythical ideal Instrument would have properties 
somewhere between AXAF and XTE, 

Consider an X-ray source with fluxes observed with a detector with 
collecting area A c but a detector area of Aj and a background flux of b. 

If the source Increases Its flux by As then In time t the excess counts 
will be AsA c t. The total counts will be sA c t+bA d t and so if one requires 
an n a detection then 


AsA c t > n/(sA c + bA d )t (E.l) 


Solving for t gives 
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where r ■ A ( j/A c . For a proportional counter, r * 1, but for an Imaging 
system A^ Is the area of one pixel and r can be quite small. A simple form 
of Equation (E.2) was used by Gorensteln (1979). 

For an Imaging system with r « 1 the Internal background, br, Is 
effectively zero. Equation (E.2) reduces to 


‘in, > 7S7— <E-3> 

45 A 1m 

where subscript 1m denotes the Imaging case. Notice that larger areas are 
preferred. Since Imaging experiments are more costly to build, they tend 
to have smaller collecting areas. The detector with the larger area will 
out perform a smaller detector, with no Internal background, for all 
sources with s > s^. It Is easy to show that 






1m 


(E.4) 


which Is Independent of as and n. Comparing XTE (br = 5 x IQ" 3 ) with AXAF, 
one finds that XTE Is more sensitive to variability for all sources with 
flux > 10“ 3 ct/cn^-sec or roughly 1 x 10“^ erg/cm 2 -sec which Includes all 
sources considered in this thesis. Note: this has assumed sources with 
similar fluxes will give similar count rates for the two detectors. This 
Is not entirely accurate since the two detectors have different band passes 
and spectral responses. 

There Is a practical minimum count rate due to the fact that an 
observatory needs to observe a large number of sources. In effect this 
means that typical observations should last a few hours. A good solid 
result requires a minimum number of photons, say 1000 (which can be divided 
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Into 10 bins with a 10% statistical error per bin). For these numbers the 
minimum count rate is ~ .1 ct/sec. Note that a source flux of 10”^ 
ct/cm 2 -sec will give 10 ct/sec in the XTE detectors. This clearly shows 
that proportional counters are "background limited" for weak sources. It 
also shows that practical considerations only require that the background 
be reduced by a factor of 100. Further reductions in the background are 
needed only to study the very low flux sources. The small number of counts 
observed from these low flux sources will be almost useless for any serious 
work on variability or spectra. 

In summary: 

1. Improved timing results requires detectors with large area. 

2. Low background Is needed to study faint sources. 

3. Instruments should be built and flown that fill the gap between XTE 
and AXAF. These instruments could be built along the lines of 
Goddard's Broad Band X-Ray Telescope (BBXRT) or the Harvard/SAO 
Large Area Modular Array of Reflectors (LAMAR). 
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APPENDIX F 

The 730 Sec Flare From NGC 4151 

Tananbaum et al. (1978) observed "significant flaring from NGC 4151 

with as much as a factor of 10 Increase (in X-ray flux) in a time as short 

as 730 sec." To date this observation has remained unconfirmed. In this 
section I reconsider the evidence for flaring and show that a smaller, 
slower change In source flux can account for the data. 

The observation was made with Uhuru which did not point at the 
source but rather scanned over It. For every pass across NGC 4151, 
Tananbaum et al. included the central 1.73 sec of data "In an attempt to 
maximize the signal from the source while minimizing the background 
counts". Although It is quite likely that different scans resulted In 
slightly different exposures, the authors give no indication as to how 
different. In addition, by effectively Ignoring all data between the 

scans, they have Ignored any changes In the detector background which will 

greatly affect their calculated probability. 

If we assume these effects are small then we can reproduce their 
numbers. For their Figure 1, they plotted 47 data bins with a mean rate of 
19.3 ct per bln of which ~ 15.7 are background. One scan had a rate of 37 
ct. Using Poisson statistics, the probability of getting 37 or more counts 
for a mean of 19.3 Is Pp ( >37 ; u=19.3) = 2.1 x 10"^ which agrees with 2 x 
1G“ 4 which they reported. However their 47 observations were made over a 
24 hour period and we now know that the flux from NGC 4151 can double In as 
little as 12 hours (Mushotzky et_ _al_. 1978, Lawrence 1980). Therefore, the 
mean for the entire observation might not be an accurate estimate for the 
flux when the flare occurred. 
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Welsskopf and Sutherland (1978) pointed out calculations, such as 
made above, are very sensitive to the estimate of the mean. To see this 
assume that the flux from NGC 4151 had doubled. The true mean Is now 
estimated to be 22.9 which corresponds to ~ 20% Increase In the total 
rate. Again using Poisson statistics Pp(>37; p*22.9) « 3.7 x 10"3. Thus a 
20% change In the mean, increases the probability that the flare was due to 
chance by a factor of 18. In addition, 3.7 x 10“3 Is the probability that 
one bln Is high. If we Include the fact that 47 bins were examined the 
probability of seeing such a deviation drops to 6%. Finally, we ask 
whether a high mean Is Inconsistent with the downward fluctuations. The 
rate 730 sec before the flare was 16. The probabilities are given by 

P {<16; ii=19.3) = .27 
P 

P (<16; u=22.9) = .084 

r 

So although the probability of getting such a downward excursion has 
dropped by a factor of 3, such an excursion Is expected for 1 bln in 12 by 
Poisson statistics alone. 

In summary, reexamination of the Uhuru observation of NGC 4151 

shows: 

1) The observation is consistent with variability on a time scale of 
less than 1 day. 

2) Slow changes In the mean can greatly increase the probability of 
seeing "flares". 

3) Small changes in the background (or even the source exposure) 
would also Increase the chances of seeing "flares". 

4) In light of the above there Is no strong evidence for 700 sec 


flares. 


150 


REFERENCES 

Bailey, T.A., Smith, A., and Turnes, 1978, Nucl. Instrum, Meth., 

155, 177. 

Baity, W.A., Mushotzky, R.F., Worrall, D.M., Rothschild, R.E., Tennant, 

A.F. , and Primlni, F.A., 1983, submitted to Ap. d. 

Baity, W.A., et al., 1981 v Ap. J., 244 , 429. 

Barr, P., White, N.E., Sanford, P.W., and Ives, J.C., 1977, M.N.R.A.S., 

181, 43p. 

Bevington, P.R., 1969. Data Reduction and Error Analysis for the Physical 
Sciences (New York: McGraw-Hill Book Company). 

Boldt. E.A., Holt, S.S., Rothschild, R., and Serlemltsos, P., 1975, Proc. 

Internat. Conf. on X-Rays in Space, _1_, 69. 

Box, G.E.P., and Jenkins, G.M., 1976, Time Series Analysis: Forcastirg and 
Control (San Francisco: Holden-Day). 

Cavaliere, A., and Morrison, P., 1980, Ap. J. (Letters), 238, L63. 

Cavallo, G., and Rees, M.J., 1978, M.N.R.A.S., 183, 359. 

Coe, M.J., Bowring, S.R., Hall, C.J., and Stephen, J.B., 1983, Astrophy. 
Let., J3_, 17. 

Cooper, B.F.C. , Price, R.M., and Cole, D.J., 1965, Austr.. J. of Phys.,_18, 
589. 

Cowle, L.L., Ostriker, J.P., and Stark, A. A., 1978, Ap. 0., 226 , 1041. 
Delvaille, J.P., Epstein, A., and Schnopper, H.W., 1978, Ap. J. (Letters), 
219 , L81. 

Dower, R.G., Griffiths, R.E, , Bradt, H.V., Doxsey, R.E., and Johnson, M.D., 
1980, Ap. J., 235, 355. 

Efanov, V.A., Moiseev, I.G., Nesterov, N.S., and Shakhovskoy, N.M., 1977, 




151 


Nature, 269, 493. 

Elvis, M., 1976, M.N.R.A.S., JJ7, 7p. 

Fabian, A.C., and Rees, M.J., 1979, In Proc. I AC/Cospar Symposium on X-ray 
Astronomy, ed. W.A. Baity and L.E. Peterson (London: Peryamon 
Press), p. 381. 

Glass, I.S., 1979, M.N.R.A.S. , _186, 29p. 

Gorensteln, P., 1979, In Compact Galactic X-ray Sources, ed. by F. Lamb and 
D. Pines, Unlv, of Illinois at Urbana-Champalgn, p. 279. 

Gorensteln, P., and Mlcklewlcz, $., 1968, Rev. Scl. Instrum,, 39, 816. 

Greensteln, J.L., and Matthews, T.A., 1963, Nature, 197 , 1041. 

Gullbert, P.W. , Ross, R.R., and Fabian, A.C., 1982, M.N.R.A.S., 199, 763. 

Halpern, J.P., 1982, Ph.D. Thesis, Harvard. 

Jenkins, G.M., and Watts, D.C., 1968, Spectral Analysis and Its 
Applications (San Francisco: Holden-Day). 

Jones, T.W., Leung, C.M., Gould, R.J., and Stein, W.A. , 1977, Ap. J., 212, 
52. 

Khachlklan, E. Ye., and Weedman, D.W., 1971, Astroflzlka, ]_, 389. 

Khachlklan, E. Ye., and Weedman, D.W., 1974, Ap. J.,^92, 581. 

Klkuchl, S., Tabara, Ho, Mikaml, Y., Kawano, N., Kawajira, N., Ojlma, T., 
Tomlnu, K., Dalshldo, T. , and Kono, M., 1973, Publ. Astron. Soc. 

Jap. , 25 , 555. 

Krolik, J.K., and London, R.A., 1983, Ap. J., 267 , 18. 

Lawrence, A., 1980, M.N.R.A.S., 192, 83. 

Lawrence, A., and Elvis, M.S., 1982, Ap. J., 256 , 410. 

Lawrence, A., Giles, A.B., McHardy, I.M. , and Cooke, B.A. , 1981, 

M.N.R.A.S., _195, 149. 

LawVence, A., Pye, P.J., and Elvis, M., 1977, M.N.R.A.S., 181 , 93p. 



152 


Letter, 0., 1980, Astr* Ap., 89, 370. 

Llghtman, A.P., 1982, Ap. J., 253, 842. 

Llghtman, A.P., Giacconl, R., and Tananbaum, H., 1978, Ap. J., 224, 375. 

Llghtman, A.P., and Ryblckl, G.B., 1979, Ap. J., 232 , 882. 

Lynden-Bell, D., 1969, Nature, 223, 690. 

Lyutyi, V.M., 1977, Astr. Zh.,_54, 1153. 

Malkan, M.A., and Sargent, W.L.W., 1982, Ap. J., 254, 22. 

Marshan, F.E., Boldt, E.A., Holt, S.S., Mushotzky, R.F., Pravdo, S.H., 
Rothschild, R.E., and Serlemltsos, P.J., 1979, Ap. J. Suppl, 40, 

607. 

Marshall, F.E., Holt, S.S., Mushotzky, R.F., and Becker, R.H., 1983, Ap. J. 
(Letters) , 269, L31. 

Marshall, F.J., and Clark, G.W., 1981, Ap. J., 245 , 840. 

Marshall, N., Warwick, R.S., and Pounds, K.A., 1981, M.N.R.A.S., 194, 987. 

Mason, I.M., and Culhane, J.L., 1983, IEEE Tran. Nuc. Scl., JO, 485. 

Matllsky, T., Shrader, C., and Tananbaum, H., 1982, Ap. J. (Letters), 258, 
LI. 

Mcllwaln, C.E., 1961, J. Geophys. Res., jft, 3681. 

Mushotzky, R.F., 1983, to appear In 11th Texas Symposium on Relat i vistic 
Astrophysics . 

Mushotzky, R.F., Holt, S.S., and Serlemltsos, P.J., 1978, Ap. J. (Letters), 
225, L115. 

Mushotj^y, R.F., Marshall, F.E., Boldt, E.A., Holt, S.S., and Serlemltsos, 
P.J., 1980, Ap. J., 235,, 361. 

Oort, (I.H., 1977, Ann. Rev. Astron. Astrophys., 15, 295. 

Osterbrock, D.E., 1978, Physlca Scrlpta, ITj 285. 


Osterbrock, D.E., 1981, In Active Galactic Nuclei, ed. C. Hazard and S. 
Mltton (Cambridge: Cambridge), pp 25-50. 


153 


Pounds, K.A., 1979, Proc. Roy. Soc. London, A, 366, 375. 

Rees, M.I., Begelman, M.C., and Blandford, R.D., 1981, In 1 0th Texas 

Symposium on Relativistic Astrophysics, ed. R. Ramaty and F.E. Jones 
(New York: New York Academy of Sciences), p. 254. 

Rleke, G.H., 1978, Ap. J., 226 , 550, 

Rothschild, R., et al., 1979, Space Scl. Instrumentation, 4_, 265. 

Sanders, R.H., 1981, Nature, _294_, 427. 

Scargle, J.D., 1981, Ap. J. Supp. Series, 45, 1. 

Schmidt, M., 1963, Nature, _197_, 1040. 

Shapiro, S.L., Llghtman, A.P., and Eardley, D.M., 1976, Ap. J., 204 , 187. 
Smart, D.F., and She-?, M.A., 1967, J.6.R. , 72, 3447. 

Stein, W.A., and Weedman, D.W., 1976, Ap. J., 205 , 44. 

Tananbaum, H., 1980, In X-Ray Astronomy, ed. R. Glacconl and G. Settl 
(Dordrecht: Reldel). 

Tananbaum, H., et al., 1979, Ap. J. (Letters), 234, L9. 

Tananbaum, H., Peters, G., Forman, W., Glacconl, R., Jones, C., and Avnl, 
Y., 1978, Ap. J., ZZ3j 74. 

Tennant, A.F., Mushotzky, R.F., Boldt, E.A., and Swank, J.H., 1981, Ap. J., 
251, 15. 

Tennant, A.F., and Mushotzky, R.F., 1983, Ap. J., 264 , 92. 

Tennant, A.F., and Swank, J.H., 1983, In preparation. 

Thorne, K.S., 1974, Ap. J.,_19U 507. 

Van Bueren, H.G., 1978, Astron. Astrophys., 70, 707. 

Weedman, D.A., 1977, In Ann. Rev, of Astron. Astrophys, pp 69-95. 


154 


Welsskopf, M.C., Kahn, S.M., and Sutherland, P.6,, 1975, Ap, J. (Letters), 
199 , L147. 

Welsskopf, M.C., and Sutherland, P.6., 1978, Ap. J., 221, 228. 

Wilson, A.S., Pooley, G.G., Willis, A.G., and Clements, E.D., 1980, Ap. J. 
(Letters), 237, L61. 

Wilson, A.S., and Willis, A.G., 1980, Ap. J., _240, 429. 

Winkler, P.E., and White, A.E., 1975, Ap. J. (Letters), 199, L139. 
Wolstencroft, R.D., Gilmore, G., and Williams, P.J., 1982, M.N.R.A.S., 201, 
479, 

Yee, H.K.C., 1980, Ap. J., 2Al> 894. 



